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SECTION  I 


INTRODUCTION 


Engine  components  have  traditionally  been  designed  using  a  crack  initiation  criterion. 
This  approach  has  been  very  successful  from  a  safety  standpoint,  but  the  conservatism  inherent 
in  this  method  has  resulted  in  poor  utilization  of  the  intrinsic  life  of  the  component.  The 
development  of  high  temperature  fracture  mechanics  has  permitted  basing  residual  life  analyses 
on  a  crack  propagation  criterion.  Using  this  approach,  lives  from  a  specifed  defect  size  can  be 
calculated  and  combined  with  periodic  inspection  to  determine  component  retirement. 

Propagation  analyses  classically  employed  in  residual  life  predictions  are  deterministically 
baaed.  They  typically  account  for  materials  scattor  by  the  incorporation  of  a  Bafety  factor.  A 
more  rigorous  treatment  of  materials  scatter  is  necessary  to  permit  maximum  utilization  of 
component  life. 

The  objective  of  this  program  wae  to  develop  a  methodology  capable  of  incorporating  the 
scatter  observed  in  crack  growth  data  into  component  residual  life  analyses.  Such  a  methodology 
is  desired  for  use  in  damage-tolerant  and  retirement-for-cause  (RFC)  concepts.  This  objective 
included:  identification  of  the  distribution  functions  which  beBt  describe  crack  growth  rate 
(da/dn)  behavior;  characterization  of  fatigue-crack  progagation  (FCP)  controlling  parameters  as 
to  their  effects  on  crack  growth  rate  variability;  examination  of  the  correlations  between  this 
variability  and  propagation  life  distributions;  and  development  of  a  generic  methodology 
applicable  to  all  engine  materials,  although  particular  program  emphasis  was  placed  on  IN100, 
Waspaloy,  and  Ti  6-2-4-6. 

Two  fr^  ire  mechanics-based  statistical  models  for  fatigue  crack  growth  damage 
accumulation  in  engine  materials  were  proposed  and  investigated  (Reference  1).  These  models 
were  based  on  hyperbolic  sine  crack  growth  rate  functions  developed  by  Pratt  &  Whitney 
Aircraft  (Reference  2).  Test  results  of  IN100  (a  superalloy  used  in  the  F100  engine)  at  various 
elevated  temperatures,  loading  frequencies,  stress  ratios,  etc.,  were  compiled  and  analyzed 
statistically.  The  statistical  distributions  of  (1)  crack  growth  rate;  (2)  propagation  Life  to  reach 
any  given  crack  size;  and  (3)  crack  size  at  any  sorviee  life,  were  derived  and  reported.  It  was 
demonstrated  that  the  correlation  between  the  IN  100  test  results  and  the  two  statistical  models 
was  very  good. 

In  this  report,  the  lognormal  statistical  model  (Reference  1)  is  extended  to  the  case  of 
spectrum  loading.  Tost  results  for  IN  100  under  block  type  loading  are  generated  and  used  to 
evaluate  the  capability  of  this  model  for  predicting  the  statistical  distributions  of  (1)  the  life  to 
reach  any  given  crack  size;  and  (2)  the  crack  size  at  any  given  service  life.  A  good  correlation  is 
shown  between  the  model  and  the  test  results. 

The  lognormal  statistical  model  is  also  applied  to  a  power  law  (linear)  crack  growth  rate 
function  (i.e.,  Paris  equation).  Test  results  for  IN100  under  various  conditions  used  to  evaluate 
the  hyperbolic  sine  function  (Reference  1)  are  also  used  to  evaluate  the  Paris  function.  This 
model  is  only  applicable  within  a  certain  range  of  stress  intensity  and  care  should  be  exercised  in 
its  application.  The  Paris  model  proves  to  be  mathematically  simple  for  practical  engineering 
applications, 

Available  fatigue  crack  growth  rate  data  for  Ti  6-2-4-6  and  WaBpaloy  under  various  test 
conditions  are  also  analyzed  statistically.  The  lognormal  statistical  model,  along  with  either  the 
Paris  or  hyperbolic  sine  crack  growth  rate  function,  is  shown  to  correlate  well  with  the  test 
results. 
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The  lognormal  model,  using  both  the  hyperbolic  sine  and  the  Paris  equations,  is  applied  to 
the  homogeneous  data  set  of  Virkler,  et  al.  (References  3  and  4).  It  is  shown  that  the  Paris 
relation  does  not  provide  a  reasonable  representation  of  this  data  set.  Finally,  a  new  statistical 
theory  is  proposed  for  the  analysis  of  fatigue  crack  propagation.  This  theory  is  based  on  the 
concept  of  fracture  mechanics  and  random  processes.  Examples  are  presented  to  demonstrate 
the  application  of  the  new  theory  using  available  test  results. 
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8ECTI0N  II 


FATIGUE  CRACK  GROWTH  UNDER  SPECTRUM  LOADING  CONDITIONS 


Two  statistical  fracture  mechanics-based  models  for  predicting  the  fatigue  crack  propaga¬ 
tion  life  of  engine  materials  under  any  single  test  condition  have  been  proposed  and  investigated 
in  Reference  1.  These  models  are  based  on  the  hyperbolic  Bine  crack  growth  rate  function 
developed  by  Pratt  &  Whitney  Aircraft.  It  has  been  demonstrated  that  the  correlations  between 
the  proposed  statistical  models  and  the  test  results  are  very  good. 

A  single  test  condition  is  defined  by  a  single  constant  value  of  each  of  the  following 
parameters:  temperature  T,  loading  frequency  v,  stress  ratio  R  and  maximum  load  P„...  Service 
loading  spectra  to  engine  components,  however,  consist  of  variable  histories  of  temperature  T, 
frequency  v,  loading  magnitude  R  and  P»„,  and  holding  time  Th.  The  objectives  of  this  section 
are:  (1)  to  extend  a  statistical  model  proposed  in  Reference  1  for  any  single  test  condition  to  be 
capable  of  predicting  the  fatigue  crack  propagation  of  engine  materials  under  spectrum  loadings, 
(2)  to  conduct  experimental  tests  using  IN100  compact  tension  specimens  subject  to  loading 
spectra  in  order  to  generate  statistically  meaningful  data,  and  (3)  to  correlate  the  test  results 
with  the  proposed  statistical  model  for  verification  purposes. 

1.  Theoretical  Model 

) 

The  lognormal  crack  growth  rate  model  proposed  in  Reference  1  for  a  single  test  condition 
is  expressed  as 


Y  -  C,  sinh[C,(X  +  C*)l  +  C,  +  Z 

in  which  C,  is  a  material  constant,  C3,  C3  and  C4  are  functions  of  test  conditions  (T,  v,  R)  and 

Y-  >««  -air  •  x-iogiK  w 

where  AK  is  the  stress  intensity  range.  In  Equation  1,  Z  is  a  normal  random  variable  with  zero 
mean  and  standard  deviation  ax  that  also  depends  on  the  test  conditions  (see  Reference  l). 

In  order  to  extend  the  statistical  model  given  by  Equation  1  to  account  for  the  crack 
propagation  under  spectrum  loadings,  the  composition  of  loading  spectra  will  be  described  as 
follows:  A  loading  spectrum  is  assume  ,  to  consist  of  repeated  identical  cycle  blocks  (or  missions 
or  duty  cycles).  Each  cycle  block  (or  duty  cycle)  is  composed  of  m  different  segments  (or  m  test 
conditions).  Each  segment  or  test  condition  is  defined  by  constant  values  of  temperature  T, 
loading  frequency  v,  stress  ratio  R,  maximum  load  Pm..  and  numbers  of  load  cycles  n.  Therefore, 
the  jth  segment  in  one  cycle  block  (or  duty  cycle)  can  be  denoted  by  (T,,  Sj,  Rj,  l»j)  whereas 
one  cycle  block  consists  of  (TJt  i>j,  R)t  Pjotl,  n()  for  j=l,2,.,.,m,  as  schematically  shown  in  Figure  1. 

For  the  jth  test  condition,  or  segment,  in  one  cycle  block,  it  follows  from  Equation  1  that 
Yj  -  C,sinh[Cai(X  +  CaJ)]  4-  C4J  +  Zj  ;  j-l,2,...,m 

(oj 

in  which  Z,  is  a  normal  random  variable  with  zero  mean  and  standard  deviation  Valuee  of  C^, 

C31,  Chanel  tr^for  different  test  conditions  are  given  in  Reference  1. 
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Cycle*,  n 
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Figure  1 .  Representation  of  One  Cycle  Block 

Since  Zj  is  a  statistical  variable,  a  specimen  having  a  faster  crack  growth  rate  will  preserve 
it  over  tho  entire  crack  propagation  life  in  the  jth  test  condition,  and  the  variability  of  Zj  is 
exclusively  contributed  by  the  material  variability  (loading  iB  deterministic)  as  described  in 

Reference  1.  Hence,  under  spectrum  loadings  Z}  for  j=l,2,3 . m  are  completely  correlated  not 

only  in  one  cycle  block  but  also  over  the  entire  crack  propagation  life.  Consequently,  Equation  5 
can  be  integrated  segment-by-segment  not  only  in  one  cycle  block  but  also  over  the  entire 
propagation  life  to  obtain  the  crack  length  a(n)  as  a  function  of  the  number  of  load  cycles,  n  (or 
cycle  blocks). 

For  simplicity,  the  crack  propagation  in  the  first  cycle  block  is  illustrated  by  the  following. 
Let  ziy  be  the  y  percentile  of  the  normal  random  variable  Zj  associated  with  the  jth  segment. 

7%  -  P[Zj>*h]  -  1  -  *(xh/<rlj) 

(4) 

in  which  *  (  )  is  the  standardized  normal  distribution  function  and  a ^  is  the  standard  deviation 
of  Zj.  The  determination  of  <r^j  was  described  in  Reference  1. 

The  7  percentile,  z^,  of  Zj  can  be  obtained  from  Equation  4  as 
zh  -  (1  -  7%)  for  j-l,2,...,m 

(5) 


where  <t>  1  ( )  is  the  inverse  standardized  normal  distribution  function. 
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The  7  percentile  of  the  crack  growth  rate  for  the  jth  segment,  denoted  by  Y^,  follows  from 
Equation  3  as 

Yh  -  C,  sinh  tOj)  (X  4-  C3,)]  +  C4,  +  zh  for  j  —  1 ,2 . m 

W 

Then,  Equation  6  can  be  integrated  segment-by-segment  to  yield  the  7  percentile  of  the  crack 
length  ar  (n)  versus  the  number  of  load  cycles  n  as  follows: 

a,  (n,)  -  a0  +  2l  Aa,,k 

k-1 


a7(nt  +  na)  -  ay  (n,)  4*  3^  Aaa,k 
ay  (2  n,  )  -aT(s‘nj)  +  S  Aa.,, 

in  which  ^  is  the  initial  crack  size,  a^  (nt)  is  the  7  percentile  crack  length  after  the  first  segment 
that  consists  of  nl  load  cycles,  and  Aal  k  is  the  increment  of  the  crack  length  during  the  kth  load 
cycle  in  the  first  segment  which  is  computed  numerically  using  Equation  6  for  j«l  as  follows: 

Ylr  -  C.sinh  [CS1(X  +  C,.)1  +  C41  +  zly 

In  Equation  8,  z,  is  determined  from  Equation  5  with  j-1.  Other  notations  appearing  in 
Equation  7  are  self-explanatory  for  segment  by  segment  numerical  integrations. 

Repeating  numerioal  integrations  similar  to  Equation  7  and  using  appropriate  initial  crack 
lengths  for  each  segment  and  cycle  block,  we  obtain  one  7  percentile  for  the  crack  length,  a  (n), 
as  a  function  of  load  cycles  n.  Furthermore,  by  varying  values  of  7,  we  obtain  the  distribution  of 
the  crack  length  as  a  function  of  load  cycles  as  shown  in  Figure  2.  It  Bhould  be  noted  that  the 
numerical  integration  to  generate  a  set  of  a  (n)  is  deterministic,  straightforward,  and  very 
simple. 


The  approach  described  above  recognizes  the  difference  of  the  statistical  dispersion  of  the 
crack  growth  rate,  da/dn,  for  each  segment  (or  test  condition).  However,  since  the  spectrum 
loading  is  idealized  by  repeated  identical  cycle  blocks,  it  may  be  reasonable  to  average  the 
statistical  dispersion  of  the  crack  growth  rate  for  each  segment  over  one  cycle  block  and  to 
approximate  the  statistical  dispersion  of  the  crack  growth  rate  for  one  cycle  block  by  such  an 
average  value.  Such  an  approximate  approach  is  described  in  the  following: 

The  coefficient  of  variation  of  the  crack  growth  rate  V.  in  the  jth  segment  is  related  to  <r , 
through 


(9) 


Then,  the  coefficient  of  variation,  V,  of  the  crock  growth  rate  for  one  cycle  block  is 
approximated  by  the  weighted  average  of  Vj  associated  with  all  segments  in  one  cycle  block,  i.e. , 


V  -  2  Vj  nj  /  2  n, 

)-i  j-i 


(10) 
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and  Equation  3  is  approximated  by 


Yj  -  C,  sinh  [C„(X  +  C„)]  +  C4i  +  Z 

in  which  Z  is  a  normal  random  variable  with  zero  mean  and  standard  deviation  at  which  is 
related  to  V  (Equation  10)  through 

<r,  “  >/E(T+TO«nlO 
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Figure  2.  y  Percentiles  of  Crack  Length  Versus  Number  of  Cycles 

Therefore,  the  statistical  dispersion  of  the  crack  growth  rate  over  the  entire  loading 
spectrum  is  characterized,  in  approximation,  by  the  random  variable  Z, 

Applying  the  same  procedures  described  previously,  with  all  Z,  (j  =•  l,2,...,m)  being  replaced 
by  Z,  one  dan  obtain  the  distribution  of  the  crack  length  as  a  function  of  load  cycles  shown  in 
Figure  2. 

To  account  for  the  scatter  of  crack  propagation  life  in  the  prediction  of  engine  components, 
two  distributions  are  of  practical  importance:  (1)  the  distribution  FNUl,(n)  of  the  propagation 
life  N(a1)  to  reach  any  given  crack  length  a4  (including  the  critical  crack  length),  and  (2)  the 
distribution,  F,(n)(x),  of  the  crack  length  a(n)  after  any  specific  number,  n,  of  load  cycles.  In  fact, 
the  distribution  of  the  crack  length  as  a  function  of  load  cycles  shown  in  Figure  2  contains  all  the 
information  mentioned  above.  For  instance,  by  drawing  a  horizontal  lino  through  any  specific 
crack  length  at  in  Figure  2,  one  obtains  the  distribution,  FN{il)(n),  of  the  number  of  load  cycles 
to  reach  that  crack  length.  Likewise,  the  distribution,  F,(nl(x),  of  the  crack  length  a(n)  after  a 
given  number,  n,  of  load  cycles  can  be  obtained  by  drawing  a  vertical  line  through  n  in  Figure  2. 

A  computer  program  based  on  the  above  theoretical  model  (Equations  6-8)  has  been 
established  to  predict  the  crack  growth  damage  accumulation  for  any  given  loading  spectrum. 
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No  attempt  was  made  in  this  program  to  account  for  retardation  or  acceleration  effects, 
although  the  statistical  model  developed  can  be  modified  to  take  these  effects  into  account. 

2.  Experimental  Program 

An  experimental  test  program  has  been  carried  out  in  order  to  verify  the  capability  of  the 
proposed  statistical  model  in  predicting  the  crack  propagation  behavior  of  engine  materials 
under  spectrum  loading.  The  specimens  used  were  ASTM  compact  tension  specimens  with  a 
thickness  (B)  of  0.5  inches  (12.7  mm)  and  width  (W)  of  2.5  incheB  (63.5  mm).  The  initial  flaw 
size  (ag)  is  0.5  incheB  (12.7  mm)  and  the  final  flaw  size  (af)  is  2.0  inch  (50.8  mm). 

The  spectrum  loading  employed  was  composed  of  repeated  block  loading  (or  cycle  block) 
consisting  of  four  test  conditions  (segments!  as  shown  in  Figure  3.  Each  segment  was  applied  for 
1,000  cycles  with  the  next  test  condition  in  soccession.  The  entire  block  was  repeated  after  test 
condition  number  4.  This  procedure  was  continued  until  specimen  fracture  occurred. 


Test 

Test 

Test 

Test 

Condition 

Condition 

Condition 

Condition 

1 

2 

3 

4 

Cycles,  n 


Test  Condition  1  1200°F,  10  cpm,  R  »  0.1  1,000  Cycles 

Test  Condition  2  1350°F,  10  cpm,  R  =  0.1  1,000  Cycles 

Test  Condition  3  1200°F,  10  cpm,  R  -  0.5  1,000  Cycles 

Test  Condition  4  1000°F,  10  cpm,  R  **  0,1  1,000  Cycles 


FD  8473*4 

Figure  3.  Laboratory  Spectrum  Test  Loading  Conditions  for  IN  100 
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In  one  block  loading,  the  maximum  load,  P^,  for  each  segment  was  fixed  at  2.4  kips 
(10.7  kn)  with  all  positive  stress  ratios,  R,  in  order  to  eliminate  the  need  for  consideration  of 
retardation  or  acceleration  effects  due  to  overloads  or  underloads.  The  results  of  twelve  tests  are 
summarized  in  Table  1,  and  displayed  in  Figure  4.  A  complete  listing  of  these  tests  is  shown  in 
the  Appendix.  The  three  solid  curves  shown  in  Figure  4  are  the  theoretical  predictions  that  will 
be  described  later. 

TABLE  1.  VERIFICATION  TEST  RESULTS  FOR 
IN100 


Specimen 

No. 

Final 

Crack 

Length 

»/ 

Cyclct  to 
Failure 

N 

Cyclee  at 
1,4  inches 

1. 

2386 

1.74 

100,000 

97,600 

2. 

2387 

1.68 

66,000 

62,000 

3. 

2398 

1.68 

68,000 

66,800 

4. 

2388 

1.43 

73,000 

72,800 

S. 

2400 

1.64 

70,000 

68,600 

6. 

2401 

1.48 

97,000 

96,600 

7. 

2402 

1.63 

86,000 

82,000 

8. 

2403 

1.63 

76,000 

73,300 

e. 

2404 

1,58 

116,000 

113,600 

10. 

2406 

1.41 

81,000 

81,000 

11. 

2406 

1.63 

81,000 

80,000 

12. 

2407 

1.46 

88,000 

84,600 

Failure: 

Maan  Life  (N)  -  82,900 

Standard  Deviation  “  14,800 

■ 

Maan  Lift  at  1.4  in.  • 
Standard  Deviation  - 

-  81,500 

16,000 

For  theoretical  predictions  the  parameter  values  Caj,  C3J  and  C4).  as  well  as  the  standard 
deviation  o^,,  of  Zj  under  each  test  condition  (or  segment)  in  one  block  loading,  Figure  3,  have 
been  estimated  using  the  method  of  maximum  likelihood  in  Reference  1.  The  results  are  shown 
in  Table  2  along  with  the  total  number  of  test  specimens  used.  The  parameter  C,  is  a  material 
constant  which  is  equal  to  0.5  for  IN1Q0,  Note  that  the  parameter  values  presented  in  Table  2 
were  estimated  UBing  compiled  IN100  crack  growth  rate  data,  which  wore  generated  over  a  long 
period  of  time  for  different  purposes  (Reference  1).  The  total  number  of  test  specimens  under 
each  teat  condition  is  very  small.  Even  with  such  a  small  sample  size,  the  specimen  dimensions, 
the  initial  and  final  flaw  sizes,  and  the  maximum  cyclic  load  for  each  test  specimen  are  all 
different.  Thus,  the  data  base  used  for  estimating  the  parameter  values  is  highly  nonhomoge- 
tieous;  however,  it  is  representative  of  the  expected  data  base  found  in  industry  (see 
Reference  1). 


Using  the  statistical  model  described  in  the  previous  section,  Equation  7,  and  Table  1,  we 
construct  variouB  y  percentiles  for  the  crack  size,  a7(n),  as  a  function  of  the  number  of  loud 
cycles  n.  The  results  are  shown  in  Figure  6  as  solid  curves  and  the  numerical  vtlues  are 
presented  in  the  Appendix.  For  instance,  the  curve  associated  with  7 -10  indicates  that  the 
probability  is  10%  that  a  specimen  will  have  a  crack  size  growing  faster  than  that  denoted  by 
the  curve.  These  solid  curves  are  the  predicted  distribution  of  crack  growth  damage 
accumulation  based  on  the  statistical  model,  The  crack  propagation  associated  with  7=6,  60  and 
95  percentiles  are  also  shown  in  Figure  4  as  solid  curves  for  the  purpose  of  comparison  with  test 
results. 


Verification  Test  Results  for  IN  100 


TABLE  2. 


PARAMETER  VALUES  C*  C3  AND  C4  AS  WELL  AS 
STANDARD  DEVIATION  <r. 


Segment 

Teit 

Condition 

c, 

c, 

c* 

Number  o/ 
Specimen* 

1 

5 

3.8982 

-1.6378 

-3.9341 

0.1028 

9 

2 

2 

4.9323 

-1.4073 

-3.9895 

0.1692 

4 

3 

3 

4.3093 

-1.3032 

-4.4460 

0.1240 

4 

4 

1 

3.8033 

-1.5239 

-4.3663 

0.1673 

6 

Figure  5.  Predicted  Crack  Size  Versus  Cycles  Behavior  Under  Spectrum  Loading 

St.  Correlation  Botwoan  Statistical  Modal  and  Experimental  Taat  Rasulta 

As  observed  from  Figure  4,  the  statistical  model  predicts  faster  average  crack  growth 
damage  accumulation  and  slightly  larger  statistical  dispersion.  Thus,  from  the  analysis  and 
design  standpoint,  the  statistical  model  is  conservative.  The  correlation  between  the  test  results 
and  theoretical  prediction  appears  to  be  very  reasonable. 

The  theoretical  predictions  for  the  statistical  distribution  of  the  number  of  load  cycles  to 
reach  crack  sizes  0.8  and  1.4  inches  (20.3  and  36,0  mm)  are  presented  as  solid  curves  in  Figures  6 
and  7,  respectively.  Also  plotted  in  these  figures  as  circles  are  the  experimental  test  results  for 
twelve  specimens.  Furthermore,  predictions  for  the  probability  of  crack  exceedance,  which  is 
equal  to  one  minus  the  distribution  function  of  the  crack  size,  at  24,000,  and  36,000  load  cycles 
are  displayed  in  Figures  8  and  9  as  solid  curves,  respectively,  Also  shown  in  these  figures  as 
circles  are  the  corresponding  experimental  test  results.  Again,  Figures  6  to  9  indicate  that  the 
statistical  model  is  slightly  conservative,  in  the  sense  that  the  model  predicts  shorter  crack 
propagation  life  and  larger  statistical  dispersion.  The  correlation  between  the  theoretical  model 
and  the  test  results  is  very  reasonable. 
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TO  246064 


Figure  6 .  Distribution  of  Cycles  to  Reach  0.8  Inch,  for  the  IN 100  Mission  Verification 
Test 


F0  244045 


Figure  7.  Distribution  of  Cycles  to  Reach  1.4  Inches  for  the  IN100  Mission  Verification 


Figure  8 


0.5  O.fl  0.7  0.8  0.0  1.0 

Crack  Siza  a  -  In. 
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Probability  of  Crack  Exceedance  at  24,000  Cycles  for  IN  100  Mission 
Verification  Testing 


Figure  9.  Probability  of  Crack  Exceedance  at  35,000  Cycles  for  1N100  Mission 
Verification  Testing 


FD  266067 


Since  the  initial  crack  size  is  deterministic,  which  is  equal  to  0,5  inch  (12.7  mm),  the 
deviation  between  the  theoretical  predictions  and  the  test  results  is  expected  to  increase  as  the 
number  of  load  cycles  increases.  This  is  clearly  demonstrated  by  the  crack  exceedance  curves 
shown  in  Figures  8  and  9.  Furthermore,  the  statistical  dispersion  of  the  crack  size  increases  with 
respect  to  the  number  of  load  cycles  as  exemplified  by  these  figures,  Hence,  the  maximum 
discrepancy  occurs  when  the  random  number  of  load  cycles  reaches  a  1.4  inch  (35.6  mm)  crack, 
that  is  very  close  to  the  critical  crack  size  (see  Figure  4).  It  is  observed  from  Figure  7  that  the 
discrepancy  in  predicting  the  average  number  of  load  cycles  to  reach  a  1.4  inch  (35.6  mm)  crack 
is  approximately  13%. 

Some  of  this  conservatism  may  be  due  to  retardation  at  the  crack  tip  caused  by  the  change 
in  yield  stress  as  a  function  of  temperature.  The  size  of  the  plastic  zone  is  inversely  proportional 
to  the  square  of  the  yield  stress,  At  the  1350°  F  test  condition,  the  lower  yield  stress  will  cause  a 
larger  plastic  zone  which  will  retard  the  crack  when  the  temperature  is  decreased  to  120O°F. 
This  will  also  occur  from  the  1200°F  to  the  1000°F  temperature  change.  The  initial  retardation 
of  the  crack  growth  rate  in  each  test  condition  would  result  in  a  longer  test  life  than  the 
predicted  life,  as  the  computer  program  (theoretical  model)  does  not  address  retardation  effects. 

From  the  statistical  standpoint,  the  main  reason  for  a  slight  discrepancy  comes  from  the 
compiled  data  base  as  described  in  the  following. 

The  theoretical  prediction  is  computed  based  on  the  parameter  values  C^,  C3,,  and  crtJ 
for  j  =  1,  2,  3,  4.  These  parameters  values  were  estimated  using  a  compiled  data  set  for  each  test 
condition  which  was  generated  over  a  long  period  of  time  for  different  purposes.  Ao  a  result,  the 
specimen  preparations,  heat  treatments  and  processes  of  materials,  test  machines,  measure¬ 
ments  and  environments,  etc.  will  not  be  identical  to  those  specimens  tested  under  the  present 
program.  The  number  of  test  specimens  compiled  for  the  data  base  for  each  test  condition  is  also 
very  small  (see  Table  1).  Likewise,  the  data  base  is  highly  inhomogeneous,  indicating  that  the 
specimen  dimension,  initial  and  final  flaw  sizes,  and  the  maximum  applied  load  are  all  different 
for  each  test  specimen.  Consequently,  it  may  be  reasonable  to  expect  Borne  inconsistencies 
between  the  highly  inhomogeneous  data  base,  from  which  the  theoretical  model  parameters  are 
calibrated,  and  the  homogeneous  test  results  generated  in  the  present  program.  Nevertheless, 
the  worst  discrepancy  is  only  13%  on  the  conservative  side,  which  is  quite  reasonable  in 
practice. 

In  practical  situations,  however,  the  data  base  is  usually  not  plentiful,  highly  nonhomoge- 
neous,  and  compiled  over  a  long  period  of  time  under  different  test  environments.  Nevertheless, 
such  a  homogeneous  and  limited  data  base  can  be  used  to  calibrate  the  model  parameters.,  and 
the  present  statistical  model  is  capable  of  predicting  the  crack  propagation  behavior  under 
spectrum  loading. 

One  important  aspect  of  the  fracture  mechanics  approach  is  that  available  crack  growth 
rate  data  from  inhomogeneous  data  bases  can  be  pooled  together  to  estimate  the  crack  growth 
rate  parameters,  although  the  type  of  specimen,  the  specimen  dimension,  initial  and  final  flaw 
lengths,  and  the  maximum  cyclic  load  for  each  test  specimen  may  be  different  in  the  data  base. 
In  contrast  to  other  statistical  models  (e.g.,  References  14,  15,  18-21),  the  statistical  model 
proposed  herein  is  based  on  fracture  mechanics,  thus  it  possesses  the  same  important  advantage 
pertinent  to  the  fracture  mechanics  approach 

Finally,  the  present  statistical  model  does  not  address  the  retardation  or  acceleration 
effects.  Further  research  is  needed  in  this  regard  in  order  to  provide  a  better  prediction 
capability  necessary  for  retirement-for-cause  analysis  of  engine  components  under  service 
loading  spectra. 


SECTION  III 


LOGNORMAL  STATISTICAL  MODEL  APPLIED  TO  THE  PARIS  EQUATION  FOR  IN100 


The  lognormal  ciack  growth  rate  model  has  been  applied  to  data  for  IN  100  at  various  test 
conditions  using  the  hyperbolic  sine  (SINH)  function  (Reference  1).  In  thiB  section,  the 
lognormal  model  is  developed  for  the  Paris  crack  growth  rate  function.  Emphasis  is  on  the 
simplicity  of  this  model  for  practical  application.  The  distributions  of  propagation  life  to  reach 
any  specific  crack  size  and  the  distribution  of  crack  size  at  any  service  life  are  derived.  The 
correlation  between  the  extrapolated  test  results  for  1N100  and  the  statistical  model  is  shown  to 
be  reasonable. 


The  Paris  crack  growth  rate  equation  is  given  by 

'  Q'AK|‘  (13) 

in  which  a  =  crack  size  after  n  load  cycles,  AK  =  stress  intensity  range,  Q  and  b  are  functions  of 
temperature  T,  loading  frequency  v,  stress  ratio  R  and  others. 

To  account  for  the  statistical  variability  of  the  crack  growth  rate  da(n)/dn,  the 
deterministic  crack  growth  rate  function  (Equation  13)  iB  randomized  as  follows 

•3T  -  X(AK)Q|4K1‘  (u) 

in  which  X(AK)  is  a  non-negative  random  function  taking  values  around  unity,  and  introduced 
to  reflect  the  statistical  nature  of  the  crack  growth  rate. 

Taking  the  logarithm  on  both  sides  of  Equation  14,  one  obtains 


Y  -  bU  +  q,  +  Z(AK) 


(IB) 


in  which 

Y  “  ,o«inr-  u  -  lo«AK 

q  -  log  Q  ,  Z(AK)  -  log  X(AK) 


(16) 


Whereas  X(AK)  is  a  random  function,  the  solution  for  the  statistical  distribution  of  the 
crack  growth  damage  accumulation  is  rather  complicated  (e.g.,  Reference  5).  In  this  connection, 
two  extreme  cases  of  the  random  function  X(AK)  should  be  mentioned.  At  one  extreme,  X(AK) 
is  completely  independent  at  any  two  values  of  the  stress  intensity  range,  referred  to  as  the 
white  noise  process.  Based  on  the  central  limit  theorem,  it  can  be  shown  that  the  statistical 
variability  of  the  crack  size  or  the  fatigue  life,  after  integrating  Equation  14,  is  the  smallest 
within  the  class  of  random  functions,  Hence,  it  is  unconservative  for  engineering  analyses  and 
applications.  Also,  the  mathematical  solution  for  the  white  noise  random  process  model  is 
difficult.  At  another  extreme,  the  random  function  X(AK)  is  totally  correlated  at  any  two  values 
of  the  stress  intensity  range,  indicating  that  X(AK)  is  a  random  variable,  i.e,,  X(AK)  =  X.  For 
the  case  of  random  variable  X,  the  statistical  dispersion  uf  the  fatigue  crack  growth  damago 
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accumulation  is  the  largest  in  the  class  of  random  functions.  Consequently,  the  random  variable 
model  is  conservative  for  the  prediction  of  the  crack  propagation  life  and  for  practical 
applications. 

Because  of  mathematical  simplicity  in  analysis  and  the  conservative  nature  in  crack 
propagation  prediction,  we  shall  investigate  the  random  variable  model,  i.e.,  X(AK)  -  X  is  a 
positive  random  variable  taking  values  around  unity.  It  is  assumed  that  X  follows  the  lognormal 
distribution  with  a  median  equal  to  unity.  Tests  for  goodness  of  fit  for  the  lognormal 
distribution  will  be  performed  later. 

It  follows  from  Equation  16  that  Z  =  log  X  is  a  normal  random  variable  with  zero  mean.  By 
virtue  of  Equation  15  the  log  crack  growth  rate,  Y  =  log(da/dn),  is  a  normal  random  variable 
with  mean  value  gyand  standard  deviation  <xy  given  by 


Hy  =>  bU  +  q 


(17) 


(18) 


in  which  it  is  obvious  that  the  mean  value  ^yis  a  function  of  AK,  and  the  standard  deviation  ayis 
a  constant. 


The  parameters  b  and  q  (or  Q),  as  well  as  the  standard  deviation  <rt  (or  <ry),  given  in 
Equations  15  through  18  car  be  estimated  from  the  test  results  of  the  crack  growth  rate  versus 
the  stress  intensity  range  using  Equation  15  and  the  method  of  maximum  likelihood  estimate. 
Since  Y  and  Z  are  normal  random  variables  and  Equation  3  is  linear,  the  method  of  maximum 
likelihood  is  identical  to  the  linear  regression  analysis  and  the  method  of  least  squares. 

Teat  results  from  compact  tension  specimens,  for  log  crack  growth  rate  Y  =  log  (da/dn) 
versus  the  log  stress  intensity  range  U  =  log  AK  are  presented  in  Figure  10  as  discrete  points  for 
five  IN100  specimens  in  the  test  condition  No,  1  (T  =  1000°F,  v  =  10  cpm,  R  =  0.1).  The  crack 
growth  rate  data  shown  in  Figure  8  have  been  used  to  estimate  b,  Q  and  <ry  using  the  method  of 
maximum  likelihood.  The  r  suits  are  presented  in  the  first  row  of  Table  3  and  plotted  in  Figure 
10  as  a  solid  straight  line  for  ay  =  0.  The  maximum  likelihood  estimates  of  b,  Q  and  <ry, 
determined  from  the  test  results  of  Y  versus  U  using  ASTM  CT  specimens,  are  shown  in  Table  3 
for  different  test  conditions  (v,R,T). 

The  power  law  for  the  crack  growth  rate  (Equation  14)  holds  only  in  the  central  region  of 
AK.  In  the  regions  where  AK  is  either  too  low  or  too  high,  the  crack  growth  rate  appears  to 
behave  asymptotically.  Hence  several  data  points  for  the  crack  growth  rate  have  already  been 
censored  in  Figure  10.  Also,  the  crack  growth  rate  data  seem  to  scutter  around  the  solid  straight 
line  indicating  the  validity  of  the  Paris  function. 

To  show  the  validity  of  the  assumption  that  Z  follows  the  normal  distribution,  sample 
values  of  Z,  denoted  by  z.,  are  computed  from  test  results  of  the  log  crack  growth  rate  Y  =  log 
da/dn  versus  log  stress  intensity  range  U  -  log  AK,  denoted  by  (yj,Uj),  using  Equation  15 

z,  =  y;  -  buj  -q  for  j  =  1 ,2,...,n 

in  which  b  and  q  have  been  estimated  by  the  method  of  maximum  likelihood  and  n  is  the  total 
number  of  test.  data. 
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Figure  10.  Data  for  Crack  Growth  Hate  Versus  Stress  Intensity  Range  for  Test 
Condition  No.  1 

TABLE  3.  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  b, 

Q,  STANDARD  DEVIATION  <ry  AND 
COEFFICIENT  OF  VARIATION  OF  da/dn 


Tut 

Coe/,  of 

No.  of 

Condition*  b 

Q 

Variation 

Data 

1  2.0617 

5.010X10”* 

0.1616 

36.0% 

261 

2  2.9216 

7.409X10-' 

0.1465 

34.5% 

147 

3  2.8439 

7.141X10-' 

0.1038 

24.3% 

129 

4  2.3860 

1.870X10"* 

0.1961 

46,2% 

185 

6  2.3180 

3.141X10”' 

0.1131 

20.5% 

338 

Average 

33.3% 

•Test  Condition! 

1.  T  -  1000*  F,  »  -  10  cpm,  R  -  0,1 

2.  T  -  1360*  F,  v  -  10  cpm,  R  -  0.1 

3.  T  -  1200*  F,  r  -  10 

cpm,  R  -  0.6 

4,  T  -  1200*F,  >  -  20 

cpm,  R  -  0,05 

8.  T  -  1200°  F,  v  -  10 

com,  R  -  0,1 

Sample  data  of  Z  (j=l,  2, ....  n)  for  test  condition  No.  3  are  plotted  on  normal  probability 
paper  in  Figure  11  as  circles  along  with  a  straight  line  representing  the  estimated  normal 
distribution  for  Z  (with  zero  mean  and  standaid  deviation  of  determined  previously,  see  Table 
3).  A  linear  scale  is  used  in  Figure  11  in  which  the  sample  data  Zj  a,e  arranged  in  an  ascending 
order,  i.e.,  z,  <  z2  :£  z3...  z„,  and  the  ordinate  corresponding  to  Zjis  given  by  <t>_l[j/(n4 1)]  with 
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4>_1  (  )  being  the  inversed  standardized  normal  distribution  function.  It  is  observed  from  Figure 
11  that  the  sample  values  of  Z  are  scattered  around  the  straight  line  without  a  nonlinear  trend, 
indicating  that  the  normal  distribution  is  acceptable. 


Figure  11.  Normal  Probability  Plot  for  Z  (Teat  Condition  No.  3) 


The  Kolmogorov-Smimov  test  for  goodness  of  fit  was  performed  to  determine  the  observed 
K-S  statistic  Dn.  The  normal  distribution  was  found  to  be  acceptable  at  least  at  a  10%  level  of 
significance. 


Since  Z  -  log  X  and  Y  -  log  da/dn  are  normal  random  variables,  the  crack  growth  rate  G  = 
da/dn  follows  the  lognormal  distribution.  The  coefficient  of  variation,  V,  of  da/dn  is  related  to 
the  standard  deviation,  <ry  =  aM,  through 


(20) 


The  coefficients  of  variation,  V,  for  the  crack  growth  rate,  G  =  da/dn,  for  five  test  conditions  are 
also  presented  in  Table  3. 


The  distribution  of  Z,  and  Y  =  log[da/dn],  are  both  normal  with  the  same  standard 
deviation,  <rt  -  ay,  that  has  been  determined  previously.  Let  zy  be  the  y  percentile  of  Z,  i.e., 


y%  -  P[Z>zv]  -  1  -  *(zy/<rt) 


(21) 


or  inversely, 

z7  -  <r,  4>-I(l  -  y%) 


(22) 
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The  7  percentile  of  the  log  crack  growth  rate  Y,  denoted  by  yy  (AK,b,q),  follows  from  Equation  3 
as 


yT(AK,b,q)  -  bU  +  q  + 

in  which  zy  is  given  by  Equation  22.  Note  that  y7  (AK,b,q)  is  a  function  of  b  and  q  which  in  turn 
depend  on  test  conditions  (r,  R,  and  T). 

Therefore,  by  varying  the  value  of  7,  one  obtains  the  distribution  of  the  log  crack  growth 
rate  in  terms  of  percentiles.  The  results  are  shown  in  Figure  12  for  test  condition  No.  1  (v  = 
10  cpm,  R  =  0.1  and  T  =  1000°F).  As  an  example,  the  crack  growth  rate  path  associated  with  7 
=  10  indicates  that  the  probability  is  10%  that  a  specimen  will  have  a  growth  rate  faster  than 
that  shown  by  the  curve. 


log(delta  k) 

FO  256034 


Figure  12.  Percentiles  of  Log  Crack  Growth  Rate  as  a  Function  of  Log  Stress  Intensity 
Range  for  Test  Condition  No.  1 


The  7  percentile  of  the  random  variable  X,  denoted  by  X  ,  is  computed  from  Equation  4  as 


X,  -  (10)** 


(24) 


in  which  z,f  is  given  by  Equation  10. 


The  7  percentile  of  the  crack  size  after  n  load  cycles,  denoted  by  a^,  is  obtained  by 
substituting  Equation  24  into  Equation  4 


da 

TiT 


XtQ1AK/ 


(25) 
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in  which  the  stress  intensity  range  AK7  for  the  ASTM  CT  specimen  is  expressed  as 


AK7 


AP 

B  y'W 


flaT] 


where  B  =  specimen  thickness,  W  =  specimen  width,  AP  =  applied  load  range  and 

f[aJ  -  -  2--“.-7r  (0.866  +  4. 64«  —  13.32aa+ 14 .72a8— 5.6a4) 

(1  —  a) ' 

“  ~w 


(26) 


(27) 

(28) 


Thus,  Equation  25  can  be  integrated  numerically  over  particular  limits  to  obtain  a  set  of 
crack  lengths,  a7(n),  versus  the  number  of  load  cycles,  n,  for  different  y  percentiles.  It  should  be 
mentioned  that  the  numerical  integration  using  Equation  25  for  each  y  percentile  of  the  crack 
length,  aY(n),  is  deterministic  and  straightforward. 

Test  environments,  including  the  initial  crack  size  the  final  crack  Bize  a,,  the  maximum 
load  P^,,  and  the  width  W  and  thickness  B  of  the  ASTM  CT  specimen,  are  assumed  for  each 
test  condition  ns  Bhown  in  Table  4.  The  necessity  to  assume  homogeneous  test  environments  in 
Table  4  for  the  correlation  study  has  been  described  previously  (Reference  1).  Integrating 
Equation  25  and  using  Tables  3  and  4,  one  obtains  five  sets  of  crack  length,  a>(  as  a  function  of 
the  number  of  load  cycles,  n,  for  different  y  percentiles.  Only  the  results  for  test  conditions 
No.  1  and  No.  2  are  displayed  in  Figure  13. 


TABLE  4.  AS3UMED  HOMOGENEOUS  TEST  ENVIRON¬ 
MENTS  FOR  TEST  SPECIMENS 


Tnt 

Condition 

(in.) 

«/ 

(in.) 

W 

(in.) 

B 

(in.) 

(hipt) 

l 

0.5 

2.0 

2.502 

0.250 

1.400 

2 

0.5 

2.0 

2.503 

0.501 

2.600 

3 

0,5 

1.8 

2.000 

0.500 

2.600 

4 

0.5 

2.0 

2.506 

0.500 

3.000 

6 

0.5 

2.0 

2.497 

0.601 

3.633 

»„  “  Initial  clack  aUc 

a,  “  final  crack  alia 

W  -  apaciman  width 

B  -  apaciman  thicknaaa 

P„,  —  maximum  load 


In  the  life  prediction  of  engine  components,  two  statistical  distributions  are  of  practical 
importance;  (1)  the  distribution  function,  FN(i  v,  of  the  number  of  cycles,  N(a,),  to  reach  any 
given  cruck  length,  and  (2)  the  distribution  ‘function,  F,(x),  of  the  crack  length  a  after  any 
number  of  load  cycles  n.  Since  Figure  13  represents  the  distribution  of  the  crack  length  as  a 
function  load  cycles  n,  it  contains  all  the  information  needed  to  determine  the  distributions 
mentioned  above.  For  instance,  by  drawing  a  horizontal  line  in  Figure  13  through  a  crack  length 
of  interest,  the  distribution  for  the  number  of  cycles  to  reach  that  crack  length  is  obtained. 
Likewise  drawing  a  vertical  line  in  Figure  13  through  a  given  load  cycle  n,  one  obtains  the 
distribution  of  the  crack  length  after  n  load  cycles.  The  complement,  F*(x),  of  the  distribution 
function,  Ft(x),  of  the  crack  length,  i.e.,  F|(x)  -  1  —  F,(x),  is  the  probability  that  the  crack 
length  after  n  cycles  will  exceed  a  certain  value  x.  Hence,  the  plot  of  F*(x)  is  referred  to  as  the 
crack  exceedance  curve. 
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Crack  Size,  a 


For  instance,  the  distribution  function,  FN(t  ),  for  the  number  of  load  cycles  to  reach  the 
crack  lengths  =  1.0  and  2.0  inches  is  obtained  from  Figure  13  and  displayed  in  Figure  14  as  a 
solid  curve  for  test  condition  No.  1.  For  the  same  condition,  the  crack  exceedance  curve  F*(x) 
after  n  =  25,000  cycles  is  shown  in  Figure  15  as  a  solid  curve. 

In  a  similar  manner,  the  distribution  functions  for  the  number  of  load  cycles  to  reach 
certain  crack  sizes  and  the  exceedance  curves  after  certain  number  of  load  cycles  for  test 
conditions  No.  2  to  No.  5  are  presented  in  Figures  16  through  23  as  solid  curves.  Note  that  the 
solid  curves  depicted  in  Figures  14  through  23  are  obtained  based  on  the  statistical  model  that 
utilizes  only  the  crack  growth  rate  data,  e.g.,  Figure  10,  for  estimating  the  model  parameters  b,  Q 
and  <ry  =  <r,  as  shown  in  Table  3. 

1.  Correlation  with  IN100  Teat  Results 

A  qualitative  correlation  Btudy  of  the  IN  100  results  is  carried  out  in  the  following  manner. 

Test  results  of  the  crack  growth  rate  for  each  specimen  are  best  fitted  by  the  hyperbolic 
sine  crack  growth  rate  function  proposed  by  Pratt  &  Whitney  Aircraft  (Reference  2). 

Y  -  C,  sinh  [CS(U  +  C,)]  +  C4 

(29) 

in  which  Y  and  U  are  given  in  Equation  16,  Cj  is  a  material  constant,  and  C2,  C3  and  C4  are 
constants  depending  on  the  test  condition  (v,  R,  T).  Equation  29  involving  four  parameters  is 
supposed  to  provide  a  better  fit  to  the  crack  growth  rate  data. 

Given  the  crack  growth  rate  data,  the  least  squares  best  fit  procedures  to  estimate  values  of 
Cj,  C3  and  C4  for  each  specimen  were  described  in  Reference  2  and  the  results  of  C2,  C3  and  C4  for 
each  specimen  in  five  test  conditions  were  presented  in  Reference  1.  Then,  Equations  29  and  16 
can  be  integrated  to  yield  the  crack  size,  a,  aa  a  function  of  load  cycles,  n.  It  has  been  shown  that 
the  crack  growth  damage  accumulation  a  thus  reproduced  by  integrating  Equations  29  and  16 
correlates  very  well  with  the  test  results  for  each  specimen  (Reference  1). 

In  order  to  correlate  the  nonhomogeneous  test  results  with  the  predictions  based  on  the 
statistical  model  proposed  previously,  homogeneous  test  environments  have  been  assumed  for 
each  test  condition,  as  shown  in  Table  4.  The  maximum  loads  Pmtl  given  in  Table  4  are  chosen  in 
order  to  avoid  excessive  extrapolation  far  into  the  region  of  AK  in  which  actual  test  results  do 
not  exist. 

By  integrating  Equations  16  and  29  over  the  homogeneous  test  environments,  with 
appropriate  values  of  C,,  C2,  C3  and  C4  for  each  specimen  (see  Reference  1),  one  obtains  five  sets 
of  the  homogeneous  crack  size,  a,  as  a  function  of  load  cycles  n.  These  homogeneous  data  sets  are 
referred  to  as  the  extrapolated  test  results,  since  they  are  not  the  results  obtained  directly  from 
experimental  tests.  The  extrapolated  test  results  for  the  crack  size  a  versus  load  cycles  n  for  the 
five  test  conditions  are  given  in  Reference  1  and  only  the  results  for  test  conditions  No.  1  and 
No.  2  are  depicted  in  Figure  24.  A  comparison  between  Figures  13  and  24  indicates  that  the 
correlation  between  the  theoretical  model  and  the  extrapolated  test  results  is  reasonable. 
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Figure  14.  Distribution  of  Cycles  To  Reach  Given  Crack  Size  for  Test  Condition  No.  I 
(a)  a  1.0  Inch,  and  (b)  a  -  2.0  Inches 
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Crack  Size,  a  -  In. 


Figure  15.  Crack  Exceedance  Curve  After  25,000  Cycles  for  Test  Condition  No.  1 


From  the  extrapolated  teat  results  for  the  crack  size  a  versus  the  number  of  load  cycles  n, 
one  obtains  (1)  extrapolated  test  data  for  the  number  of  load  cycles  to  reach  any  given  crack  size 
a,  by  drawing  a  horizontal  line  through  a1(  and  (2)  extrapolated  test  data  for  the  crack  size  after 
any  number  of  load  cycles  nt  by  drawing  a  vertical  line  through  nt. 

The  distribution  function  is  constructed  from  the  extrapolated  test  data  obtained  above  by 
arranging  them  in  an  ascending  order.  The  ordinate  of  the  ith  data  is  given  by  i/(m+l)  where  m 
is  the  total  number  of  data  points.  For  instance,  one  obtains  five  data  points  from  Figure  22  by 
drawing  a  horizontal  line  through  at  ■*  1.0  inch.  These  five  data  points,  for  the  number  of  cycles 
to  reach  a  crack  size  of  a,  *=  1  inch,  are  ranked  in  an  ascending  order  where  the  ordinate  of  the 
ith  data  point  is  equal  to  i/6  (i=l,  2, ...,  5).  These  results  are  shown  in  Figure  14(a)  by  circles.  In 
like  manner,  the  distributions  of  the  extrapolated  test  results  for  the  number  of  load  cycles  to 
reach  other  crack  sizes,  as  well  as  their  crack  exceedances  after  some  numbers  of  cycles  have 
been  constructed  and  shown  in  Figures  14  through  23  as  circles.  It  is  observed  from  Figures  14 
through  23  that  the  correlation  between  the  statistical  model  (solid  curves)  and  the  extrapolated 
test  results  (circles)  is  very  reasonable. 
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Figure  17.  Crack  Exceedance  Curve  after  10,000  Cycles  for  Teat  Condition  No.  2 
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Figure  19,  Crack  Exceedance  Curve  after  40,000  Cycles  for  Test  Condition  No.  3 
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Figure  22.  Distribution  of  Cycles  to  Reach  (liven  Crack  Size  for  Test  Condition  No.  5; 
(a  I  a  -  0.8  Inch,  and  (b)  a  =  2.0  Inches 
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SECTION  IV 


ANALY8IS  OF  Tl  6-2-4-B  AND  WASPALOY  DATA 


Available  Ti  6-2-4-6  and  Waspaloy  crack  propagation  test  results  from  constant  amplitude 
cyclic  loading  with  various  test  conditions  were  compiled.  The  specimen  dimensions,  initial 
crack  size,  final  crack  size,  and  maximum  load  for  each  specimen  in  each  test  condition  are 
shown  in  Tablea  5  and  6.  The  crack  growth  rate  data  for  all  specimens  in  each  test  condition  are 
also  pooled  together  for  analysis  purposes.  The  pooled  data  for  Waspaloy  and  Ti  6-2-4-6  from 
test  condition  Number  1  are  shown  therein,  respectively. 


TABLE  5.  SPECIMEN  GEOMETRY  AND  MAXIMUM  LOAD 
FOR  EACH  TEST  SPECIMEN  OF  Ti  6-2-4-6 


Teit 

Condition 

Specimen 

No, 

“o 

(in,) 

A 

(ini 

H 

W 

(in.) 

K3SI 

WtJPHl 

'  -  Room  Temperature 

1536 

0.8138 

1.7724 

0.249 

2.496 

0.5460 

-  20  Hi 

1580 

0.8130 

1.6130 

0.2970 

2.6007 

0.7200 

1  =  0.1 

1582 

0.6880 

2,2080 

0.3090 

2.0086 

1.0640 

1662 

0.7360 

1,0620 

0.4990 

1.4970 

0.6800 

1663 

0,7580 

1,0590 

0.4990 

1.4990 

0.6810 

1676 

0.39B0 

1.1780 

0.3520 

1.9910 

0.9580 

1677 

0.4040 

1.1900 

0.3630 

1.9987 

0.9880 

1680 

0.4730 

1.1164 

0.3735 

1.9965 

1.0060 

1B72 

0.1919 

0,6025 

0.4990 

0.B999 

1.1960 

1535 

0.7835 

1.7618 

0.2485 

2.4999 

0.6090 

“  Room  Temperature 

0423 

0.3850 

0,7973 

0.0979 

2.0000 

1.890 

“  10  cpm 

0448 

0.8628 

1.7005 

0.2520 

2.5001 

0.538 

-  0.1 

1484 

1.3789 

1.7677 

0.5000 

2.4951 

1.037 

1886 

0.6786 

1.1716 

0.5004 

1.9929 

1.182 

1886 

0.6733 

1.1711 

0.6005 

2.0066 

1.138 

1892 

0.5149 

1.4487 

0.3890 

2.4993 

1.296 

2103 

0.5577 

1.2686 

0.3700 

2.5003 

1.301 

“  Rocm  Temperature 

1744 

0.0828 

0.2261 

0.2935 

0.996 

6.611 

“  10  cpm 

1741 

0.0766 

0.1982 

0.2900 

0.996 

7.157 

-  0.0 

1742 

0.0887 

0,2637 

0,2970 

0.996 

6.816 

1532 

0.115 

0,3120 

0.2950 

0.997 

4.273 

-  800-F 

0358 

0.597 

1,7926 

0.2750 

2.2520 

0.624 

-  10  opm 

0359 

0.6697 

1,5878 

0.2750 

2.2556 

0.626 

-  0.1 

0400 

0.3696 

0.8369 

0.1250 

2.0000 

4.844 

0410 

0.3674 

0.0907 

0,1205 

2.0000 

4.891 

1245 

0.5453 

0.7775 

0.5010 

1.0030 

0.306 

1881 

0.3068 

1,4964 

0,6240 

1.9959 

1.703 

1883 

0.3679 

1,4014 

0,5040 

2,0046 

1,577 

-  400° ¥ 

0413 

0.8325 

2,0600 

0.2500 

2.5000 

0.3330 

-  20  Hi 

0415 

0.3741 

0,8276 

0.1260 

2.0000 

2.4500 

-  0.1 

0416 

0.8426 

2,0489 

0.2500 

2.5002 

0.2690 

1678 

0.4260 

0,9300 

0,3440 

1.9972 

1.1160 

1886 

0.6128 

1,4764 

0,6016 

1.9986 

1.1210 
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TABLE  6.  SPECIMEN  GEOMETRY  AND  MAXIMUM  LOAD 
FOR  EACH  TEST  SPECIMEN  OF  WASPALOY 


Test 

Condition 


Specimen 

No. 


T  -  1200' F 
u  -  20  Hi 
R  -  0.06 


0736  1.0706  1.7666  0.432  2.502 

1464  0.6000  2.0820  0.500  2.500 

1470  0.8510  1.7460  0.491  2.491 

1304  1.0830  1.6670  0.500  2.509 


2 


T  -  1200'F 

1004 

0.0703 

0.3326 

0.301 

0.998 

11.513 

v  -  10  cpm 

1301 

1.1072 

1.7094 

0.601 

2.519 

1.769 

R  -  0.05 

1013 

0.9688 

1.8653 

0.132 

2.484 

1.261 

3 

T  -  800“F 

1002 

0.8339 

1.7598 

0.499 

2,605 

2.574 

u  -  10  cpm 

1003 

0.0635 

0.3047 

0.303 

0.997 

13.316 

R  -  0.06 

1018 

1.0134 

1.8944 

0.298 

2.60C 

1.417 

Using  the  hyperbolic  sine  crack  growth  rate  function  (Equation  29),  and  the  lognormal 
statistical  model  (Reference  1),  we  obtain  the  maximum  likelihood  estimates  for  C3>  C3,  C4  and 
the  standard  deviation  <ry  =  <r,.  These  results  are  presented  in  Tables  7  and  8.  The  crack  growth 
rateB  based  on  the  values  of  C3,  Ca,  C3  and  C4  given  in  these  tables,  for  Z  *  0,  are  shown  as  solid 
curves  in  Figures  25  and  26  (test  condition  Number  1). 

TABLE  7.  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  C4,  C3,  C4, 
STANDARD  DEVIATION  a  «  a,  AND  COEFFICIENT 
OF  VARIATION,  V,  OF  da/dn  FOR  Ti  6-2-4-6  (HYPER¬ 
BOLIC  SINE  FUNCTION) 


Test 

Condition 

c, 

c , 

c, 

V 

No.  of 
Data 
Points 

1 

0.7 

4.68 

-1.0898 

-5.4878 

0.0672 

16.6% 

269 

2 

0.7 

4.5715 

-1.1088 

-6,378 

0.1156 

27.1% 

179 

3 

0.7 

3.0949 

-0.9651 

-6.7946 

0.1174 

27.5% 

50 

4 

0.7 

2.1964 

-0,9086 

-5.49b8 

0.0661 

13,0% 

251 

5 

0.7 

3.6301 

-1.0641 

-6.7081 

0.0939 

21.9% 

153 

TABLE  8.  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  C2,  Ca, 
C4,  STANDARD  DEVIATION  -  <r,  AND  COEFFI¬ 
CIENT  OF  VARIATION,  V,  OF  da/dn  FOR 
WASPALOY  (HYPERBOLIC  SINE  FUNCTION) 


Test 

Condition 

c, 

e. 

c, 

V". 

V 

No.  of 

Data 

Points 

1 

0.5 

3,5484 

-1.4064 

-5.0269 

0,0783 

18.2% 

130 

2 

0.5 

4.0154 

-1,3694 

-4.8774 

0.0470 

10.9% 

60 

3 

0.6 

4.3817 

-1.4941 

-4.9847 

0.0654 

12.8% 

67 
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Crack  growth  rate  data  for  each  specimen  are  beat-fitted  by  the  hyperbolic  sine  function 
(Equation  3),  to  determine  the  parameter  values  Ca,  C3  and  C4.  The  results  of  the  parameter 
values  for  each  specimen  are  shown  in  Tables  9  and  10. 


TABLE  9.  PARAMETERS  Ca,  C3,  AND  C4  FOR  EACH  TEST 
SPECIMEN  OF  TITANIUM  WITH  C,  =  0.7,  HYPER¬ 
BOLIC  SINE  CRACK  GROWTH  RATE  FUNCTION 


Tut 

Condition 

Spteimtn  No, 

c, 

e. 

c, 

R 

Type  of 
Spteimtn 

1  — 

1636 

4.6692 

-1.0206 

-5.735 

0.1 

MOT 

1530 

4.1696 

-1.2061 

-5.1217 

0.1 

MCT 

1682 

6.0229 

-1.3336 

-4.8700 

0.1 

MCT 

1662 

4.4731 

-1.0217 

-5.6861 

0.1 

MCT 

1663 

4.7667 

-1.0399 

-6.6614 

0.1 

MCT 

1676 

4.1696 

-0.9866 

-6,8264 

0.1 

MCT 

1677 

4.2736 

-0.9868 

-5.8294 

0.1 

MCT 

1680 

4.1694 

-1.2464 

-4,9161 

0.1 

MCT 

1072 

3.7766 

-1.0662 

-5,6379 

0.1 

FRT 

1636 

4.8336 

-1,0530 

-6,6762 

0.1 

MCT 

2  — 

0423 

4.6117 

-1.1744 

-5.2811 

0.1 

ON 

0448 

4.4617 

-1.0333 

-6,6818 

0.1 

MCT 

1484 

4.8831 

-1.4262 

-4.2868 

0.1 

MCT 

1886 

6.2786 

-1.1319 

-6.2332 

0.1 

MCT 

1888 

4.0642 

-0.8937 

-6.1046 

0.1 

MCT 

1892 

4.6282 

-1.0629 

-5.6214 

0.1 

MCT 

2103 

6.1136 

-0.9827 

-6.7010 

0.1 

MCT 

3  - 

1744 

6.4666 

-1.3166 

-4.6476 

0.0 

CN 

1732 

3.1613 

-0.8446 

-6.6793 

0.0 

CN 

1741 

6.1208 

-1.2970 

-4.7369 

0.0 

CN 

1642 

4.7786 

-1.4152 

-4,3071 

0.0 

CN 

4  - 

0668 

2.1664 

-0,8647 

-6.6210 

0.1 

MCT 

0369 

1.6211 

-0.4380 

-8.1880 

0.1 

MCT 

0409 

3.0048 

-1.7464 

-3.7829 

0.1 

CN 

0410 

3.3867 

-1.4631 

-4.4611 

0.1 

CN 

1246 

2.2046 

-0.8064 

-6,6619 

0.1 

MCT 

1881 

2.7011 

-1.0981 

-6.1862 

0.1 

MCT 

1883 

2.3800 

-0.9634 

-5.4186 

0.1 

MCT 

a  — 

0413 

3.6878 

-1.1146 

-5,6747 

0.1 

MCT 

0416 

3.4328 

-1.2919 

-6.1373 

0.1 

CN 

0416 

3.8776 

-0.8233 

-6,4866 

0.1 

MCT 

1678 

3.7288 

-1.2226 

-5.2324 

0.1 

MCI’ 

1886 

3,6905 

-1.1537 

-5,3067 

0.1 

MCT 

TABLE  10. 

PARAMETERS  C.( 

C„  AND 

C,  FOR  EACH  TEST 

SPECIMEN  OF  WASPALOY  WITH  C, 

-  0.5,  HYPER- 

BOLIC  SINE  CRACK  GROWTH  RATE  FUNCTION 

Tut 

Condition 

Spteimtn  No, 

* 

C, 

R 

Typo  of 
Spteimtn 

i 

0736 

3.9744 

-1.6278 

-4.8472 

0,06 

MCT 

1464 

3.4B68 

-1,4725 

-4.8366 

0.06 

MCT 

1470 

4.4467 

-1.2349 

-6.3546 

0.05 

MCT 

1304 

3,8747 

-1.6438 

-4,7752 

00.1 

MCT 

2 

1004 

4,9009 

-1.4710 

-4.6876 

0.06 

CN 

1301 

3.7554 

-1.2698 

-5.1172 

0,05 

MCT 

1013 

3.8140 

-1.4935 

-4.6702 

0.1 

MCT 

3 

1002 

4.4664 

-1,6408 

-4.8281 

0.06 

MCT 

1003 

4.9438 

-1.5226 

-4,9012 

0.06 

CN 

J.018 

3.827 

-1.3911 

-5.2284 

0.06 

MCT 

In  order  to  perform  correlation  studies,  homogeneous  test  environments  should  be 
assumed,  since  each  teat  specimen  has  different  geometry,  initial  crack  size,  final  crack  size,  and 
maximum  load  P^,  as  shown  in  Tables  6  and  6.  Homogeneous  test  environments  assumed  for 
Waspnloy  are  shown  in  Table  11. 

TABLE  11.  ASSUMED  HOMOGENEOUS  TEST  ENVIRON¬ 
MENTS  FOR  EACH  TEST  CONDITION  FOR 
WASPALOY 


Tmt 

Condition 

»o 

_ (in.) 

<v 

(In.) 

m 

m 

mm 

mMm 

Type  of 
Specimen 

1 

2.0 

0.5 

2.5 

3.6 

MCT 

2 

0.5 

2.5 

6.2 

MCT 

3 

0.6 

0.5 

2.5 

4.0 

MCT 

The  distribution  of  the  crack  growth  damage  accumulation,  a,  as  a  function  of  load  cycles, 
n,  based  on  the  lognormal  statistical  model  can  be  obtained  using  Tables  S  and  11.  Then,  the 
distribution  functions  for  (1)  the  number  of  load  cycles  to  reach  any  specific  crack  size,  and  (2) 
the  crack  exceedance  probability  ut  any  given  number  of  load  cycles,  can  easily  be  determined 
(Reference  1 ).  The  distribution  functions  for  the  number  of  load  cycles  to  reach  crack  Bizes  0.875 
and  2.0  inches  are  displayed  as  solid  curves  in  Figures  27(a),  28(a),  and  29(a),  for  the  three  test 
conditions.  Likewise,  the  probabilities  of  crack  exceedance  at  n  -  50,000,  10,000  and  60,000 
cycles,  for  the  three  test  conditions,  are  shown  aB  solid  curves  in  Figures  27(b),  28(b),  and  29(b). 

With  the  aid  of  the  best-fitted  parameter  values  for  Ca,  C3,  and  C4,  given  in  Table  10,  the 
crack  growth  damage  accumulation,  a,  for  each  individual  specimen,  is  obtained  by  integrating 
the  crack  growth  rate  equation  and  using  the  assumed  homogeneous  teat  environments  given  in 
Tahle  11.  The  distribution  of  the  number  of  load  cycles  to  reach  any  specific  crack  size  and  the 
crack  exceedance  probability  at  any  given  number  of  load  cycles  is  then  established.  These 
extrapolated  results  are  presented  as  circles  and  triangles  in  Figures  27  through  29.  The 
homogeneous  test  environment  in  Table  11,.  is  assumed  to  avoid  excessive  extrapolation  into  the 
region  of  AK  for  which  test  data  do  not  exist.  This  may  be  difficult  to  achieve  because  the  crack 
growth  rate  data  for  each  specimen  do  not  cover  the  same  region  of  AK.  Therefore,  the 
extrapolation  cannot  be  avoided  for  some  specimens.  This  problem  becomes  more  serious  for  the 
Ti  fi-2-4-6  data  and  will  be  discussed  later. 

It  is  observed  from  Figures  27  through  29  that  the  correlation  between  the  lognormal 
statistical  model  and  the  extrapolated  tost  results  is  very  reasonable  for  WaBpaloy. 

In  a  similar  manner,  homogeneous  test  environments  are  assumed  for  Ti  6-2-4-0.  Based  on 
the  lognormal  statistical  model,  the  distribution  functions  for  the  number  of  cycles  to  reach  two 
different  crack  sizes  are  shown  as  solid  curves  in  Figures  30(a)  through  34(a)  while  the  crack 
exceedance  curves  are  shown  in  Figures  30(b)  through  34(b).  The  corresponding  extrapolated 
test  results  are  shown  in  the  same  figures  by  circles  and  triangles.  It  is  observed  that  except  for 
test  condition  number  6,  the  correlation  is  Ibbs  than  satisfactory. 
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Figure  27,  Distribution  Curve  and  Crack  Exceedance  Curve  for  Waspaloy  for  Teat 

Condition  No.  1;  (a)  Distribution  of  Cycles  to  Reach  Given  Crack  Size,  and 
(b)  Crack  Exceedance  Curve  After  50,000  Cycles 
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Figure  28.  Distribution  Curve  and  Crack  Exceedance  Curve  for  Waspaloy  for  Test 
Condition  No.  2;  (a)  Distribution  of  Cycles  to  Reach  a  Given  Crack  Size, 
and  (b)  Crack  Exceedance  Curve  After  10,000  Cycles 
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Figure  30.  Distribution  Curve  and  Crack  Exceedance  Curve  for  Titanium  for  Test 
Condition  No.  1,  Using  the  Hyperbolic  Sine  Model;  (a)  Distribution  of 
Cycles  to  Reach  Given  Crack  Size,  and  ( b )  Crack  Exceedance  Curve  After 
30,000  Cycles 
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Figure  31,  Distribution  Curve  and  Croc’  Exceedance  Curve  for  Titanium  for  Test 
Condition  No.  2,  Using  the  Hyperbolic  Sine  Model;  (a)  Distribution  of 
Cycles  to  Reach  Given  Crack  Size,  and  (b)  Crack  Exceedance  Curve  After 
20,000  Cycles 
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Figure  33,  Distribution  of  Cycles  to  Roach  Given  Crack  Size  for  Titanium  for  Test 
Condition  No,  4,  Using  the  Hyperbolic  Sine  Model;  (a)  Distribution  of 
Cycles  to  Reach  Given  Crack  Size,  and  (b)  Crack  Exceedance  Curve  After 
5,000  Cycles 
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Figure  34.  Distribution  Curve  and  Crack  Exceedance  Curve  for  Titanium  for  Test 
Condition  No.  5,  Using  the  Hyperbolic  Sine  Model;  (a)  Distribution  of 
Cycles  to  Reach  Oiven  Crack  Size,  and  (b)  Crack  Exceedance  Curve  After 
100,000  Cycles 
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A  careful  examination  of  the  data  indicates  that  the  ext.apolation  for  some  specimens  is 
excessive,  due  to  the  following  reasons:  (1)  The  crack  growth  rate  data  for  each  specimen  cover 
different  ranges  of  AK,  and  therefore,  it  is  impossible  to  assume  a  homogeneous  test 
environment  without  excessive  extrapolations  for  some  specimens,  and  (2)  The  hyperbolic  sine 
function  has  two  asymptotes,  and  the  extrapolation  is  very  sensitive  in  the  region  of  AK  where 
an  asymptote  occurs.  In  test  condition  Number  1,  for  instance,  the  crack  growth  rate  data  of  one 
specimen  (Number  7AN1592)  covers  only  a  very  small  range  of  AK.  Consequently,  the  best- 
fitted  hyperbolic  sine  function  results  when  the  left  asymptote  falls  into  the  region  of 
extrapolation  as  shown  in  Figure  35.  Therefore,  the  propagation  life  is  too  long  as  indicated  in 
Figure  30  by  numerical  values.  Consequently,  it  seems  reasonable  to  censor  this  specimen.  With 
such  a  specimen  being  censored,  the  results  are  shown  in  Figures  36(a)  and  (b).  It  is  observed 
that  the  correlation  is  very  reasonable.  Similar  situations  exist  for  test  conditions  Numbers  2 
through  4. 

The  Paris  crack  growth  model  was  also  used  to  analyze  statistically  the  Waspaloy  and 
Titanium  data. 

From  the  pooled  crack  growth  rate  data  for  each  teat  condition,  the  method  of  maximum 
likelihood  was  employed  to  estimate  the  parameters  b  and  log  Q  as  well  bb  the  standard 
deviation  <ry  =  <rE.  The  results  are  presented  in  Tables  13  and  14,  respectively,  for  Waspaloy  and 
Titanium.  Furthermore,  the  crack  growth  rate  data  for  each  test  specimen  were  best-fitted  by 
the  Paris  function  to  determine  the  corresponding  b  and  log  Q  values.  The  results  for  each 
specimen  are  shown  in  Tables  15  and  16,  respectively. 

Following  the  Bame  procedures  described  in  Section  III  and  using  the  assumed  homoge¬ 
neous  test  environments  given  in  Tables  11  and  12,  the  correlation  between  the  statistical  model 
and  the  extrapolated  test  results  was  obtained.  The  correlation  results  for  the  distribution  of 
random  number  of  load  cycles  to  reach  any  given  crack  size  and  the  crack  exceedance  curve  are 
depicted  in  Figures  37  through  44.  It  is  observed  that  the  correlation  is  very  good  for  Waspaloy 
data  sets  whereas  a  considerable  improvement  in  correlation  has  been  achieved  for  Titanium 
data  sets. 


TABLE  12.  ASSUMED  HOMOGENEOUS  TEST  ENVIRON¬ 
MENTS  FOR  EACH  TEST  CONDITION  FOR 
TITANIUM 


r««t 

Condition 

(in.) 

<•/ 

(in.) 

£ 

(in.) 

W 

(in.)  . 

Pm.x 

(kip*) 

Type  of 
Speoimtn 

1 

0.4 

1.7 

0.3830 

1.9087 

1.40 

MCT 

2 

0.4 

1.7 

0.3890 

2.4943 

2.40 

MCT 

a 

0.08 

0.4 

0.3000 

1.0000 

7.90 

CN 

4 

0.4 

1.7 

0.5000 

2.0000 

4.64 

MCT 

6 

0.4 

1.7 

0.3440 

1.9972 

1.12 

MCT 
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TABLE  13.  MAXIMUM  LIKELIHOOD 
ESTIMATE  OF  b,  LOG  Q 
AND  STANDARD  DEVIA¬ 
TION  <r,  =  <ry  FOR  WASPA- 
LOY  (PARIS  FUNCTION) 


Tut  Condition 

b 

log  Q 

»y“  »« 

1 

2,1841 

-8.0979 

0,0826 

2 

2.1429 

-7.8101 

0.0476 

3 

2.4193 

-8.5940 

0.0576 

TABLE  14.  MAXIMUM  LIKELIHOOD 
ESTIMATE  OF  b,  LOG  Q 
AND  STANDARD  DEVIA¬ 
TION  <r,  -  <rv  FOR  TITANI¬ 
UM  (PARIS  FUNCTION) 


Toil  Condition 

b 

log  Q 

a,  -  Cy 

1 

3.6062 

-9.4166 

0,0797 

2 

3.6240 

-9,2847 

0,1180 

a 

2,8691 

-8,4737 

0.1273 

4 

2.1210 

-7,5039 

0.0926 

5 

3.1094 

-9,0305 

0,0767 

TABLE  16.  PARAMETERS  b  AND  LOG  Q 
FOR  EACH  TEST  SPECIMEN  OF 
WASPALOY,  (PARIS  FUNCTION) 


Specimin  No, 

Tut  Condition  t 

b 

log  Q 

0736 

2.2398 

-8,2642 

1484 

2.2937 

-8,2108 

1470 

2.5171 

-8.4611 

1304 

2.0887 

-7,9979 

Tut  Condition  3 

Specimin  No, 

b 

log  q 

1001 

2.5790 

-8.4792 

1301 

2.1521 

-7.8828 

1013 

2.1302 

-7.7509 

Tut  Condition  3 

Specimin  No, 

5 

log  q 

1002 

2.4769 

-8.6585 

1003 

2.8654 

-9,2381 

1018 

2.1978 

-8,2906 
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TABLE  16.  PARAMETERS  b  AND  LOG  Q 
FOR  EACH  TEST  SPECIMEN  OF 
TITANIUM  (PARIS  FUNCTION) 


S pic  i nun  No, 

Teit  Condition  1 

b 

lot  Q 

1635 

3.8178 

-9.7068 

1630 

3.6888 

-9.4006 

1680 

3.6838 

-9.4074 

1682 

4.0286 

-10.0166 

1662 

3.6677 

-8.4288 

1663 

3.8253 

-9.6413 

1676 

3.6012 

-9.2797 

1677 

3.6373 

-9.4230 

1680 

3.7600 

-9.5362 

1072 

2.8632 

Test  Condition  3 

-8.7004 

Specimen  No. 

b 

lot  Q 

“* 

0423 

3.3768 

-9.2429 

0448 

3.3188 

-9,1072 

1484 

5.0674 

-11.2912 

1888 

4.0281 

-9.7903 

1886 

3.4606 

-9.2101 

1882 

3.6220 

-8.4686 

2103 

4.3748 

Test  Condition  3 

-10.0228 

Specimen  No, 

6 

lot  Q 

1744 

4.5104 

-10.6678 

1632 

3.0872 

-8.6326 

1741 

3.8822 

-8.9043 

1742 

4.8882 

Teel  Condition  4 

-11.0684 

Specimen  No, 

b 

log  q 

' 

0368 

1.0848 

-7.2854 

0368 

1.0619 

-7.3831 

0409 

2.7940 

-8.6696 

0410 

2.4115 

-7.9669 

1246 

1.7779 

-7.1221 

1831 

2.3706 

-7,8062 

1883 

2.23111 

Teet  Condition  6 

-7.6221 

Specimen  No, 

b 

lot  Q 

0413 

3.2663 

-9.2172 

0416 

2.6365 

-8.6358 

0418 

2.9630 

-8.B137 

1878 

2.9365 

-8.8684 

1883 

2.7603 

-8.6468 
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Figure  38.  Distribution  Curve  and  Crack  Exceedance  Curve  for  Wanpuloy  for  Test 

Condition  No.  2,  Using  the  Paris  Model;  (a)  Distribution  of  Cycles  to  Reach 
Given  Crack  Size,  and  ( b )  Crack  Exceedance  Curve  After  10,000  Cycles 
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Figure  41  Distribution  Curve  and  Crack  Exceedance  Curve  for  Titanium  for  Test 

Condition  No.  2,  Using  the  Paris  Model;  (a)  Distribution  of  Cycles  to  Reach 
Given  Crack  Size,  and  ( b )  Crack  Exceedance  Curve  After  20,000  Cycles 
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Figure  12. 
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Distribution  of  Cycles  to  Reach  Given  Crack  Size  for  Titanium  for  Test 
Condition  No.  3,  Using  the  Paris  Model;  (v)  Distribution  of  Cycles  to  Reach 
Given  Crack  Size,  and  ( b )  Crack  Exceedance  Curve  After  5,000  Cycles 
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Figure  43.  Distribution  of  Cycles  to  Reach  Given  Crack  Size  for  Titanium  for  Test 

Condition  No.  4,  Using  the  Paris  Model;  (a)  Distribution  of  Cycles  to  Reach 
Given  Crack  Size,  and  (b)  Crack  Exceedance  Curve  After  5,000  Cycles 
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SECTION  V 


STATISTICAL  ANALYSIS  OF  A  HOMOGENEOUS  DATA  SET 

A  homogeneous  crack  propagation  data  set  for  2024-T3  aluminum  center-cracked 
specimens  under  constant  amplitude  cyclic  loading  generated  in  References  3  and  4  is  analyzed. 
The  lognormal  statistical  model,  developed  under  the  present  program,  and  an  advanced  model 
currently  being  developed  under  another  program,  have  been  uBed.  The  results  of  such  analyses 
reveal  the  importance  of  both  the  selection  of  the  crack  growth  rate  function  and  the  data 
analysis  technique  used  for  deriving  the  crack  growth  rate  data  from  the  experimental 
measurements. 

Sixty-four  sample  functions  for  the  crack  growth  damage  accumulation,  a,  versus  the 
number  of  load  cycles,  n,  are  shown  in  Figure  46.  The  initial  half-crack  length  is  9  mm  and  the 
final  half-crack  length  is  49.8  mm.  The  crack  growth  rate  data,  da/dn,  versus  the  stress  intensity 
range,  AK,  are  derived  from  Figure  45  using  the  seven-point  polynomial  method  and  the  secant 
method.  Based  on  the  hyperbolic  sine  crack  growth  rate  function,  (Equation  1),  the  method  of 
maximum  likelihood  is  employed  to  estimate  the  parameters  C2,  C3,  C4  and  the  standard 
deviation  *  =  <r,  (see  Reference  1).  The  results  are  shown  in  Table  17.  It  is  observed  from  Table 
17  that  the  statistical  dispersion  of  the  crack  growth  rate,  da/dn,  is  smaller  when  the  seven-point 
polynomial  method  is  applied. 


Figure  45.  Homogeneous  Data  Set  of  Virkler,  et  al.  with  =  23.35  kn  and  R  =  0.2 

Using  the  lognormal  statistical  model  and  integrating  the  crack  growth  rate  equation 
(Equation  1),  for  different  y  percentiles,  one  obtains  the  distribution  of  the  half-crack  length  as 
a  function  of  load  cycles  n.  The  results  nm  shown  in  Figure  46,  in  which  the  seven-point 
polynomial  method  haB  been  used,  Note  mat  the  stress  intensity  range  AK  for  the  CCT 
specimen  is  gi  ven  by 

AK  =  Bee  (ira/2W) 
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TABLE  17.  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  C2,  C„  C4, 

STANDARD  DEVIATION  -  <r,  AND  COEFFICIENT  OF 
VARIATION,  V,  OF  da/dn;  Cj-  0.5 


Analyiit 
Ttchniqu * 

c, 

c, 

c4 

", 

V 

7-Point 

Polynomial 

3.4477 

-1.3002 

-4.5348 

0.08235 

19.13% 

Secant 

method 

5.0148 

-1.1336 

-5,0738 

0.08589 

22.36% 

Figure  46.  Lognor'mal  Statistical  Model  Prediction  of  Virkler,  et  al.  Data  Set 


A  comparison  between  Figures  45  and  46  indicates  that  the  lognormal  statistical  model 
correlates  well  with  the  tost  results  for  the  50%  crack  propagation  life.  However,  the  statistical 
disperaion  of  the  crack  growth  damage  accumulation  based  on  the  model  is  larger  than  the 
experimental  test  results. 

The  distribution  functions  for  the  number  of  cycles  to  reach  half-crack  lengths  of  21  mm 
and  49.8  mm  are  presented  in  Figures  47  and  48,  respectively.  In  these  figures,  the  solid  curve 
and  the  dashed  curve  represent  the  results  baaed  on  the  lognormal  statistical  model.  The  former 
is  obtained  using  the  7-point  polynomial  method,  whereas  the  latter  is  obtained  using  the  secant 
method.  The  experimental  test  results,  Figure  45,  are  shown  by  circles  in  Figures  47  and  48. 
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Figure  47.  Distribution  Function  of  Number  of  Cycles  to  Reach  21mm  Half-Crack 
Length 

The  probabilities  of  crack  exceedance  at  n  =  150,000  cycles  are  shown  in  Figure  49  as  a 
solid  curve,  a  dashed  curve,  and  circles  which  have  the  same  meaning  as  those  in  Figures  47  and 
48. 

As  observed  from  Figures  47  through  49,  the  correlation  between  the  lognormal  statistical 
model  and  the  test  results  is  very  good  for  the  50%  crack  propagation  life.  However,  the  model 
predicts  a  larger  statistical  dispersion  for  crack  growth  damage  accumulation  than  the  actual 
test  data  which  results  in  a  conservative  prediction  for  the  crack  propagation  life.  Furthermore, 
the  seven-point  polynomial  method  appears  to  be  superior  to  the  oocant  method. 

The  observation  that  the  lognormal  statistical  model  results  in  a  larger  statistical 
dispersion  for  the  crack  growth  damage  accumulation  was  expected  as  described  in  Reference  1. 
This  stems  from  the  assumption  that  the  crack  growth  rate  is  completely  correlated  and  hence  Z 
is  a  random  variable.  At  another  extreme,  when  the  crack  growth  rate  is  assumed  to  be 
completely  uncorrelated  at  any  two  different  values  of  log  AK,  one  obtains  a  white  noise  process. 
The  white  noise  process  model  results  in  the  smallest  statistical  dispersion  for  the  crack 
propagation  life  (Reference  1)  as  evidenced  by  the  results  presented  in  References  3  and  4. 
However,  for  the  white  noise  model,  the  prediction  for  the  crack  propagation  life  is 
unconservative. 
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Figure  48,  Distribution  Function  of  Number  of  Cycles  to  Reach  49.8mm  Half-Crack 
Length 


The  actual  test  results  in  Figure  45  indicate  a  definite  correlation  for  the  crack  growth  rate 
at  different  values  of  log  AK.  Therefore,  a  statistical  model  taking  into  account  such  a 
correlation  has  been  explored  in  Section  VI,  which  iB  referred  to  as  the  advanced  statistical 
model.  A  complete  development  for  the  analytical  solution  of  such  a  model  has  been 
accomplished  under  the  sponsorship  of  another  program  (Reference  28).  It  cun  be  shown  that 
the  statistical  dispersion  of  the  crack  propagation  life  decreases  as  the  correlation  for  the  crack 
growth  rate  reduces. 


Based  on  the  advanced  model  with  the  correlation  parameter  A  =  40,000  (see  Section  VI), 
the  distribution  functions  for  the  number  of  cycles  to  reach  half-crack  lengths  21  mm  and 
49.8  mm  are  indicated  by  dotted  and  dash-dot  curves  in  Figures  47  and  48.  The  dotted  curve 
represents  the  results  of  the  advanced  model  using  the  seven-point  polynomial  method,  while 
the  dash-dot  curve  denotes  the  results  using  the  secant  method.  Similarly,  the  crack  exceedance 
curve  at  n  *>  150,000  cycles  is  Bhown  in  Figure  49  by  dotted  curve  and  dash-dot  curve  having  the 
same  meaning  as  those  in  Figures  47  and  48. 

As  indicated  from  Figures  47  through  49,  the  advanced  model  with  the  method  of  Beven- 
point  polynomial  correlates  well  with  the  experimental  test  data  (circles).  For  the  advanced 
model,  the  7-point  polynomial  method  is  definitely  superior  to  the  secant  method. 
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Figure  49 ,  Crack  Exceedance  Curves  at  150,000  Cycles 


In  conclusion,  the  Paris  crack  growth  rate  function,  Equation  13,  has  been  used  in 
conjunction  with  both  the  lognormal  and  advanced  statistical  models.  The  correlation  between 
the  test  data  and  the  prediction  baaed  on  the  Paris  crack  growth  rate  funtion  is  extremely  poor 
even  for  the  50%  crack  propagation  life.  Therefore,  it  is  not  worthwhile  to  present  the  results. 

The  results  of  analyses  for  the  homogeneous  data  set  presented  above  indicate  two  salient 
features  associated  with  the  prediction  of  fatigue  crack  propagation.  For  a  particular  type  of 
specimen  under  a  particular  type  of  fatigue  loading,  one  should  select  a  crack  growth  rate 
function  which  can  best  describe  the  crack  growth  rate  behavior.  For  the  homogeneous  data  set 
analyzed,  the  hyperbolic  sine  function  fit  the  crack  growth  rate  data  very  well,  whereas  the  Paris 
function  did  not.  Using  the  method  of  maximum  likelihood,  the  goodness-of-fit  can  be  judged 
from  the  standard  deviation  <ry  =  A  crack  growth  rate  function  has  a  better  fit  if  the 
corresponding  standard  deviation  at  is  smaller. 

Data  for  the  crack  size,  a,  versus  the  number  of  load  cyci  a,  n,  are  measured  directly  from 
experiments.  Then,  the  crack  growth  rate  data  are  derived  from  a  versus  n  measurements. 
Various  analysis  techniques  for  obtaining  the  crack  growth  rate  data  have  been  propofod  in  the 
literature,  such  as  the  secant  method,  the  method  of  seven-point  polynomial,  etc.  Unfortunately, 
the  accuracy  for  the  statistical  prediction  of  fatigue  crack  propagation  depends  on  the  data 
reduction  methodology  employed.  Further  research  is  needed  to  identify  a  best  technique  for 
the  reduction  of  the  crack  growth  rate  data  as  well  as  the  measurements  of  a  versus  n. 
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SECTION  VI 


ADVANCED  STATISTICAL  MODEL  FOR  FATIGUE  CRACK  GROWTH 


Several  mathematical  models  have  been  proposed  for  the  prediction  of  crack  growth 
damago  accumulation  for  structures  under  dynamic  loads  based  on  the  principles  of  fracture 
mechanics  (References  6  and  7),  These  models  have  the  general  form  of 


-|f-  -  Q(K,AK,s,a,R) 

(31) 

where  a(t)  =  crack  size  at  time  t,  Q  -  a  non-negative  function,  K  =  stress  intensity  factor,  AK  ™ 
stress  intensity  range,  s  -  stress  amplitude,  R  =  stress  ratio.  For  instance,  the  well-known  Paris 
model  is  given  by 

4r  -  w4K>‘ 

(32) 

AK  ->  a(a)  As  \Zira 

(33) 

where  As  =  stress  range,  a(  )  =  a  function  depending  on  the  specimen  and  crack  geometries. 
However,  even  in  a  well-controlled  laboratory  environment,  results  obtained  from  crack  growth 
experiments  under  either  a  constant-amplitude  cyclic  loading  or  a  given  spectrum  loading 
usually  exhibit  considerable  statistical  variability.  This  iB  illustrated  in  Figure  60,  which  shows 
the  actual  crack  growth  records  of  some  fastener  hole  specimens  subjected  to  the  excitation  of  a 
specific  load  spectrum  in  a  laboratory.  Ic  is,  therefore,  not  surprising  to  find  that  statistical 
analyses  have  been  appliod  quite  frequently  to  such  problems  in  recent  years  (References  6 
through  22). 


If  we  restrict  our  attention  to  a  laboratory  setting  so  that  the  loading  time-variation  is 
deterministic,  then  the  mathematical  model,  Equation  (1),  can  be  “randomized”  as  follows: 

rtf? 

UT  “  Q(K,AK,s,a,R)  X(t) 

where  the  additional  factor  X(t)  is  a  non-negative  random  process,  and  a(t)  is  a  random  process 
representing  the  crack  length  at  time  t.  It  is  of  interest  to  note  that  Virkler,  Hillberry  and  Goel 
(References  3  and  4)  have  undertaken  simulation  studies  of  crack  propagation  which  amount  to 
assuming  X(t)  in  Equation  (34)  to  be  totally  uncorrelated  at  any  two  different  times.  At  the 
other  extreme,  Yang  (References  1,  22,  17)  has  replaced  the  random  process  X(t)  in  Equation 
(34)  with  a  random  variable,  which  is  equivalent  to  the  case  where  X(t)  is  totally  correlated  at  all 
times.  It  was  pointed  out  in  Reference  1  that  a  totally  uncorrelated  X(t)  would  lead  to  the 
smallest  statistical  dispersion  and  a  totally  correlated  X(t)  to  the  greatest  statistical  dispersion 
for  the  time  at  which  a  given  crack  size  is  reached.  A  more  realistic  modeling  of  fatigue  crack 
growth  should  lie  somewhere  between  the  two  extremes.  Therefore,  the  ability  to  account  for  an 
arbitrary  correlation  in  X(t)  is  a  major  consideration  in  the  present  section. 
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Figure  50.  Actual  Crack  Propagation  Time-Histories  of  Some  WPB  Fastener  Holes 

Prom  the  standpoint  of  fatigue  design  and  schedule  maintenance  (References  11, 14, 15), 
two  are  of  internat:  the  dretribotion  of  the  random  .rack  »M  at 

anv  riven  time  and  the  distribution  of  the  random  time  to  reach  a  given  crack  size. _  As  shown  l 

Z.»n«  i  two  problems  me  lalarrelated.  Attention  will  be  foctmad  on  the  atatotiMd 

nroDerties  of  the  latter  problem.  In  particular,  analytical  solutions  will  be  given  for  the 

statistical  moments  of  the  random  time  to  reach  any  crack  « LoMHf be 

size  A  Drocedure  to  estimate  the  parameters  in  a  power-law  mathematical  model  also  win  De 
;ZlTusing  the  fractographical  results  (Reference  23)  of  -me  74  5-T7361  aluminum 
fastener  hole  specimens  subjected  to  the  excitation  of  a  bomber  load  spectrum. 

1.  Modal  tor  X(t) 

We  shall  model  X(t)  as  a  random  pulse  train  (see  Reference  24),  i.e., 

X(t)  -  Zk  w(t,rk)  (35) 

k-1 

where  N(t)  -  a  homogeneous  Poisson  counting  process,  giving  the  total  number  of  pulses  that 
S,"  Sn  the  time* interval  (-00,  t);  r,  -  th.  arrivel  time  nf  the  kth  puleei  Z,  -  the  random 
amplitude  of  the  kth  pulse,  and 

0  <  t-r  <  A 

re  )  ’  (36) 

w(t,T)  =  w(t-r)  -  j  otherwise 


We  further  assume  that  Zk  for  different  k  are  independent,  identically  distributed  random 
variables,  with  a  common  probability  distribution  Z. 


A  typical  sample  function  of  X(t)  is  shown  in  Figure  61. 
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Figure  51.  Typical  Sample  Function  of  Random  Process  X(t) 


The  statistical  properties  of  X(t)  can  be  described  by  its  cumulant  (or  semi-invariant) 
functions.  The  mth  cumulant  function  is  given  by: 

KjX(t,),...,X(tJl  -  E(Z“]X/"'",‘""‘Ww(tI-  r)..,w(tm-  r)dr 

-  w  (37) 

in  which  E[  ]  denotes  an  ensemble  average,  X  is  the  average  arrival  rate  of  the  PoisBon  process, 
and  min(  )  indicates  the  smallest  of  the  parenthesized  quantities.  In  particular,  the  first 
cumulant  is  the  mean  function,  and  the  second  cumulant  is  the  covariance  function.  These  are 
found  to  be: 

p  -  E(X(t))  -  E[Z]X  f'  w(t  -  r)dr 

-  00 


-  E(Z]X  /"  w(u)du  -  E(Z]X  A 


(38) 
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and 


Cov[X(t,)>  X(ta)]  -  B[Zs|X/‘‘  w(t,  -  T)w(t2-  r)dr 

—  00 


*»  E(Za]X  w(u)w(t2  -  t,  +  u)du 

fe/8(l  ~  |t2- tJ/A),  |t,  -  t,|  <  A 

M>,  |t,-tJasA 


(39) 


in  which  0  =  ^E[Z2]XA. 

2.  Approximation  of  a  (t)  by  a  Markov  Random  Procaao 

We  now  re-write  Equation  (84)  as  follows! 

/]  o* 

■ar  “  Q®  if*  +  Y(t)i  (40) 

where  the  dependence  of  Q  on  K,  AK,  b  and  R  has  been  suppressed  for  simplicity.  Clearly,  Y(t) 
is  a  random  process  with  a  zero  mean  and  the  correlation  function  of  Y(t)  is  the  same  as  the 
covariance  function  of  X(t);  namely, 


(  2/3(1-  M/A),  M  <  A 

Ryy(t)  -  E(Y(t)  Y(t+r)]  -  l  (41) 

( 0,  otherwise 

A  sketch  of  this  correlation  function  is  shown  in  Figure  52.  If  the  correlation  time  of  Y  (t)  is  short 
compared  with  the  characteristic  time  of  a  (t),  then  a  (t)  is  close  to  a  diffusive  Markov  process 
(Reference  25)  which  is  governed  by  an  Ito’s  stochastic  differential  equation  (Reference  26): 

da  —  m(a,t)dt  +  <r(a,t)dB(t) 

(42) 

where  m  is  called  the  drift  coefficient,  <r  the  diffusion  coefficient,  and  B(t)  is  a  unit  Brownian 
motion  process  (also  called  the  Wiener’s  process),  which  has  the  property  that  dB(t,)  and  dB(t2) 
are  independent  for  tj  ¥>  ta. 


Tho  correlation  time  of  Y(t)  may  be  defined  us  follows: 

_  /:T|RYV(r)ldr 

Substitution  of  Equation  (41)  into  Equation  (43)  results  in  “  A/3. 


(43) 


68 


FD  256038 


Figure  52.  Autocorrelation  Function  of  Random  Process  Y(t) 


Strictly  speaking  the  a  (t)  process  In  Equation  (42)  is  an  approximation  of  that  in 
Equation  (40),  and  they  could  be  represented  more  clearly  by  two  different  symbols;  however, 
the  aame  symbol  will  be  used  in  this  report  for  both  processes  as  long  as  no  confusion  will  result. 

Stratonovich  (Reference  20)  hus  given  the  required  formulas  to  compute  the  drift  and 
diffusion  coefficients  from  the  original  physical  equation  when  the  Markov  approximation  is 
justified.  In  the  case  of  Equation  (40) 

m  -  Qm  +  /"  Q-~j- E|Y(t)Y(t  +  r)|dr 

-  Qm  +  -  q(m  +  /iA-f-) 

a  u  (44) 

<r'J  “  2  /" ..  E|Y(t)Y(t  +  r)|dr  -  2Q20A 

where  r0  >  A.  Theso  equations  imply  that  Q  and  ilQ/flu  vary  slowly  within  the  integration 
interval  to  justify  their  being  taken  outside  the  integrals.  ThuB,  the  integrals  account  basically 
for  the  contribution  towards  the  drift  and  diffusion  due  to  the  correlation  between  the  past  and 
the  present  "excitations"  Y(t).  This  contribution  Is  lumped  at  the  present,  when  m  and  a  are 
used  in  Equation  (42).  The  replacement  of  Equation  (40)  by  Equation  (42)  amounts  to 
substituting  Y(t)  by  a  white  noise.  Theoretically,  the  substitution  introduces  a  small  error 
associated  with  a  small  probability  for  da  to  become  negative.  This  error  is  negligible  so  long  as 
the  drift  dominates  the  diffusion,  which  should  be  verified  in  each  practical  case.  Stratonovich’s 
formulas  are  applicable  to  other  types  of  correlation  functions  for  Y(t)  as  long  as  Q  and  dQ/da 
vary  slowly  within  an  interval  of  r  where  such  a  correlation  is  not  negligibly  small.  In  this  case, 
the  lower  limit  of  integration  can  even  be  extended  to  ~oo.  Stratonovich’s  method  is  known  as 
the  stochastic  averaging  method,  originully  proposed  on  u  physical  ground,  but  later  proved 
rigorously  by  Khasminskii  (Reference  27)  and  given  a  rigorous  mathematical  interpretation. 
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The  transition  probability  density  qn  (a,  t  I  a0,  t0)  of  the  Markov  process  a  (t>  is  u 
conditional  probability  density  which  describes  the  distribution  of  a  (t)  under  the  condition  that 
the  initial  crack  size  is  a  (t0)  =  a0  at  an  earlier  time  t0.  It  iB  governed  by  the  following  Fokker- 
Planck  equation  (Reference  26): 


+  —  [q(m  +  P A-f-)  qj  -  -4-  (QadAq.)  -  0 


(46) 


or  by  the  adjoint  of  Equation  (46): 

x  + [«0. +#*£)]• 


«  0 


(47) 


subject  to  the  condition, 


q;(a,t0 1  Mo)  “  <Ka-a0)  (4g) 

In  Equation  (46)  Q  is  treated  as  a  function  of  "u"  which  is  a  sample  value  of  the  crack  size  a  (t), 
whereas  in  Equation  (47)  it  is  treated  as  a  function  of  a0.  These  equations  are  also  known  as  the 
Kolmogorov’s  forward  and  backward  equations,  respectively. 

3.  Time  to  Roach  a  Qivon  Crack  Sizo 

We  now  focus  our  attention  on  the  random  time  when  the  crack  size  a  (t)  reaches  a  specific 
value  a0.  This  event  may  be  considered  as  the  passage  for  the  first  time  across  an  absorbing 
boundary,  Introduce  g(a,  1 1  u0,  t0)  such  that 

g(a,  t|  a„.  t„)du  “  Prob  |u  <  a(t)  <  u  +  du  |  u0  £5  3(r)  <  aL„  t„  :£  r  <  t|, 


a,,  ^  a  <  a„ 


(49) 


Although  g  1h  similur  to  q,ai  not  a  probability  density  since  its  integration  on  a  from  a0  to  ac 
is  generally  smaller  than  one.  In  fact,  g  describes  only  the  sample  functions  which  do  not  reach 
the  absorbing  boundary  ac  before  time  t.  However,  g  also  satisfies  the  Kolmogorov  backward 
equation,  Equation  (47). 


The  integration  of  g  yields 
G(t)  -  /l'  g(a,  1 1  a„,  t(()  da 


Prob  survival  in  (t„,t)  |  u(t„)  =  a() 


(50) 


Here  the  dependence  of  G  on  Up  and  t0  has  been  suppressed.  The  G(t)  function  also  satisfies  the 
Kolmogorov  backward  equation,  numely, 


-r-  +  q(m  +  /3A-^-)  ~  +  QadA  — 4  G 

t()  Q()  «o  a„ 


0 


(51) 
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subject  to  the  following  conditions: 


G(t0)  *>  1 


(52) 


G(oo)  “  0 


(53) 


Condition  (63)  indicates  that  sooner  or  later  fatigue  failure  will  occur.  Letting  r  -  t  -  t0,  we 
obtain 


qO-.  +  -  0 


(54) 


Now  G(t0)  -  G(t)  -  1  -  G(t)  ia  the  probability  that  first  passage  (reaching  the  crack  size  ac) 
occurs  prior  to  t.  This  is  the  distribution  function  of  the  random  first  passage  time  T;  i.e., 


F v( r )  “  1  -  G(t0  +  r) 

The  probability  density  of  T  follows  from  a  differentiation  of  Equation  (66); 
Pr(r)  “  ~  3— G(t0  +  r) 


(55) 


(56) 


The  average  first  passage  time  is 

E[TJ  -  -  37-dr  -  /:G(t0+r)dr 


In  obtaining  the  second  part  of  Equation  (67)  use  has  been  made  of  the  condition  (63). 


(67) 


Integrating  Equation  (54)  and  using  Equations  (56)  and  (67),  we  obtain  am  equation  for  the 
average  first  passage  time: 


1  +  q(m  +  )-^E[T]  +  ql/?A— rE[T]  -  0 


(58) 


The  above  second  order  equation  requires  two  boundary  conditions,  one  of  which  is  clear'y 
E|T,  -  0,  if..  -  «.  m 


We  shall  assume  that  another  condition  is 


Sir E(T]  “  0,  if  a„ 


(60) 


Equations  for  higher  order  moments  can  also  be  derived  from  Equation  (54).  Note  that 


E(T“]  -  -  Sy^-dr  -  nf)  ~  '  G(t„  +  r)dr 


(61) 
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Multiplying  Equation  (64)  by  r“and  integrating  on  r,  we  obtain 

(n  +  1)E[T“1  +  Q  (  m  +  ) -a^-ElT*  + ‘]  +  Q2/3A-^-E[T“  + ‘]  -  0 

subject  to  the  conditions: 

E[Tn  +  1]-0,  ifa0-at 


a  E  T»  +  i]  »  o  if  a0  -  0 

dan  (64) 

It  is  clear  that  Equation  (64)  which  inch'd es  Equation  (60)  as  a  special  case  implies  that  a,)  *»  0  is 
a  reflective  boundary.  Equation  (62)  can  be  used  recursively  to  obtain  higher  moments  from  the 
lower  order  moments.  Equation  (62)  reduces  to  Equation  (58)  when  n  =  0. 


Note  that  Equation  (62)  is  first  order  in  dE[TB+’]/da0  which  can  be  solved  readily.  Let  Q  - 
Qof(ft0)  where  Qq  has  the  unit  of  length/time  and  f(  • )  is  a  dimensionless  function.  We  obtain 


where 


J-  -  /rEITVlexp^/^)] 


in  which  tho  common  practice  of  using  the  same  symbol  for  the  integration  variable  and  limit  of 
integration  has  been  adopted  for  convenience.  This  solution  satisfies  condition  (64).  Equation 
(66)  implies  that  d/da,,  E[T"+1]  is  nonpositive  or  E[T"+1]  is  non-increasing  when  the  initial  crack 
size  a0is  increased,  which  is  physically  reasonable. 

The  (n+l)th  momont  of  T  follows  from  integration  of  Equation  (65): 

BIT-* 'I  -  BIT- *■]<!■•  m 

which  satisfies  condition  (63). 


4.  Power-Law  Crack  Propagation 

We  now  consider  a  special  form  for  the  Q  function  which  has  been  proposed  by  several 
authors  (References  13  through  16): 


Q  -  Q0  (ff/ac)b 
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and  in  the  case  of  backward  equations,  we  have  Q  =  Q0(a0/ae)b.  In  this  case,  the  integration  in 
Equations  (66)  end  (67)  can  be  carried  out  in  closed  form.  Specifically 


*(■££).  s-i 


f  dao 


ur 


£n  u  , 


b  -  1 


(69) 


where  u  -  a0/ac.  To  compute  E[T],  we  require  J0.  From  Equation  (66)  and  with  n  =  0, 
l^-[exp  (T=F»''k)  “  H(l-b)],  bfkl 


Jo  " 


where  k  ! 


H(l-b) 


1  k 
WU> 


and  H  is  the  Heaviside  function,  i.e., 


b  -  1 


(70) 


I,  b  <  1 

0,  b  >  1 


(71) 


Equations  (69)  and  (70)  are  substituted  into  the  expression  for  dE[T]/day  which  is  then 
integrated  to  obtain  E[T];  the  results  are 


QoM 


E[TJ 


■{TzT  U-  u"l'>  +  H(l-b)  (-f)  [«ip(  -  1^5') 
”p(  -  T^T”'  '’)]}.  t>  1 


(72) 


--gir8""' 

To  compute  E(Ta]  we  require 


b  -  1 


Jl  ”  7g7/?ElTld[exp(-fr  /  -ij1)] 


(73) 
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Substituting  Equation  (72)  into  Equation  (73),  we  obtain 

Q&2  ,  r  1  1  (  k  ,-b  ,\1  (  k  ,~h\ 

“  hnr  ~  t  vt=fu  -i;Je*plT^Fu  ' 

+  H(l-b)j-J-exp(  -  T^¥-)[exp(T^-u,"b)  -l] 


-  y^-(u'“b  +  1)  -  b  ^  1 

”  uk  (-5 - finu),  b  -  1 

In  obtaining  Equation  (75),  we  have  used  the  formula 
/; U*fin u du  -  -  /;u*du/;-^- 


-  (a  +  1)  ~  8x*  +  ‘[(a  +  l)£n  x-1] 


The  second  part  of  Equation  (76)  follows  from  changing  the  order  of  integration. 

Substitute  Equations  (74)  and  (76),  respectively,  into  the  expression  for  dEfT^/da,,  and 
integrate- 


EIT-!  -  +  -E-)*1-”'-1)  -  - 


2(1  —  b)a 


(1  -uatl  -  b>) 


+  H«-b)  ji-  +  ,«p(- 

-  e*p(-  xV  “'■')[  TCXT,!  < 1  +  “ '  " l)  +  + (t^t)1|  . 


b  ¥-  1 


E[T8]  -  «nu(-5-Snu  -  -£-)<  b  “  1 


In  obtaining  Equation  (78)  we  have  used  the  formula: 

/xirSnudu“  -  Jx du /I dv 


(fin  x)8 


The  second  part  of  Equation  (79)  follows  from  an  interchange  of  the  order  of  integration. 

9.  Estimation  of  Modal  Parameters 

The  estimation  of  the  parameters  of  a  mathematical  model  will  be  illustrated  by  use  of  two 
examples.  Refer  again  to  Figure  50  which  shows  the  crack  propagation  time-histories  of  some 


1 


specimens  selected  from  the  fractographical  data  of  7476-T7351  aluminum  fastener  holes 
available  in  Reference  23.  The  data  set  is  identified  as  WPB,  indicating  that  the  specimens  were 
drilled  with  Winslow  Spacematric  machines  (W),  with  a  proper  drilling  technique  (P),  and 
subjected  to  a  given  B-l  bomber  load  spectrum  (B).  The  data  set  has  been  censored  to  include 
only  those  specimens  having  fatigue  crack  growth  through  the  crack  length  interval  from  0.004 
inch  to  0.04  inch  (corner  crack).  This  censoring  procedure  is  necessary  in  order  to  ndrmalize  the 
data  set  to  zero  life  at  0.004  inch,  thus  obtaining  a  homogeneous  crack  growth  data  base.  The 
resulting  data  set  consists  of  16  specimens.  Within  a  small  crack  size  range,  the  power-law 
propagation,  Equation  (63),  was  shown  to  be  valid  (References  13  through  17). 

The  parameters  of  our  mathematical  model,  specialized  to  a  power-law  crack  propagation 
that  fits  the  behavior  of  the  above  WPB  specimens,  have  been  estimated  using  the  experimental 
data.  We  began  by  taking  logarithms  on  the  two  sides  of  Equation  (68),  with  Q  *  Q0(a/ac)b  to 
yield 


-bloga(t)  +  log(Q0/ab)  +  log  X(t) 

(80) 

Variations  of  log  d  a  /dt  versus  log  a  were  plotted  for  the  16  specimens  in  the  WPB  data  set  as 
shown  in  Figure  53.  Implicit  in  the  mathematical  model  is  the  assumption  that  all  points  would 
fall  along  a  single  straight  line  if  the  random  element  log  X(t)  were  not  present.  Since  X(t)  has 
been  assumed  to  be  a  stationary  random  process,  the  mean  and  the  standard  deviation  of  log 
X(t)  are  constants.  A  linear  regression  analysis  was  then  carried  out  to  estimate  the  slope,  the 
intercept,  and  the  standard  deviation  of  the  random  element  log  X(t).  The  slope  is  equal  to  b, 
and  the  intercept  is  log  (Q0/  ab) . 

The  results  obtained  from  the  16  specimens  are 

b  “  0.92971, 

Q0/ab  «  1.1061  X  10“*, 

ffl«*x(t) =  0.087636 


The  linear  regression  analysis  implies  that  each  point  in  Figure  63  has  been  treated  as  an 
independent  sampling  of  a  Gaussian  random  variable.  Thus,  estimates  of  the  mean  and  variance 
of  X(t)  itself  can  be  computed  from  those  of  log  X(t),  using  the  log-normal  to  normal  conversion 
formulas;  namely, 


MX  -  exp  (  -5-  [«ri^x«)2n  lof  )  -  m 
*x  -  Mx|exp  [ffb.xwBn  10]a  -  l}  -  20 


(81) 

(82) 


Application  of  Equations  (81)  and  (82)  resulted  in 


M  -  1.0206,  0  -  0.021643 
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Figure  53.  Regression  Analysis  for  the  Estimation  of  Model  Parameters  for  WPB 
Fastener  Holes 

The  rationale  behind  these  conversion  rules  is  the  use  of  linear  regression  to  obtain  the  mean 
and  standard  deviation  for  log  X(t).  It  does  not  necessarily  require  that  X(t)  must  be  log¬ 
normal.  In  fact,  when  the  crack  size  S(t)  is  treated  as  a  diffusive  Markov  process  in  an 
approximate  sense,  the  “excitation"  process  X(t)  is  also  effectively  replaced  by  a  constant  mean 
H  plus  a  Gaussian  white  noise.  Of  course,  this  replacement  excitation  process  cannot  be  log¬ 
normal. 

Having  the  values  of  b,  Q,,,  n  and  0,  the  mean  and  mean-square  values  of  the  time  T  to 
reach  a  given  crack  size  a0can  be  computed  using  Equations  (72)  and  (78)  for  different  values  of 
k  -  aju/^AQj),  or  equivalently  for  different  values  of  A.  Figure  54  shows  the  results  of  such 
computations,  in  terms  of  E[T]  and  E(T]  ±  <rT,  where  <rT  -  (E[T!]  -  E(T)  a)1/ais  the  standard 
deviation  of  T.  The  mean  values  E[T]  are  practically  uninfluenced  by  the  choice  of  A  since  the 
last  two  terms  in  Equation  (72)  are  several  orders  of  magnitude  smaller  than  the  first  term. 
However,  the  statistical  dispersion  <rT  increases  as  A  increases.  This  A-<rT  relationship  agrees 
with  our  earlier  observation  concerning  the  results  of  "Vang  (References  1,  22),  and  of  Virkler, 
Hillberry  and  Goel  (References  3  and  4),  the  former  case  being  equivalent  to  A  — ►  oo  and  the 
latter  case  to  A  — »  0.  Estimates  of  the  mean  E[T]  and  standard  deviation,  <rT,  of  T  were  also 
obtained  directly  from  the  data  of  the  16  specimens,  and  the  results  plotted  in  Figure  55. 
Comparison  between  Figures  54  and  55  shows  that  A  *  8,000  is  a  reasonable  choice  for  this 
particular  mathematical  model  and  data  set. 
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'igure  54,  Theoretical  Mean  and  Standard  Deviation  of  Random  Time  to  Reach 

Various  Crack  Sizes  Computed  for  Different  A  Values  for  WPB  Fastener 
Holes 


Figure  55,  Mean  and  Standard  Deviation  of  Random  Time  to  Reach  Various  Crack 

Sizes  Computed  Directly  from  Soma  Actual  Time-Histories  of  WPB  Fastener 
Holes 
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Some  computer  simulated  sample  functions  using  the  present  mathematical  model 
described  in  Equations  (34)  and  (35)  are  shown  in  Figure  56.  These  have  been  obtained  using 
parameters  -  1.0206,  vxs  =  0.043286,  b  =  0.92971,  Q0/acb  =  1.1051  X  10'4,  A  =  8,000  and  an 
average  Poisson  pulse  rate  of  X  =  0.1.  Their  general  characteristics  are  remarkably  similar  to  the 
actual  records  in  Figure  50. 


Figure  56.  Simulated  Sample  Functions  of  Crack  Propagation  Time-Histories  for  WPB 
Fastener  Holes 


It  has  been  suggested  that  the  distribution  of  the  time  T  to  reach  a  given  crack  size  ac  can 
be  approximated  by  a  two-parameter  Weibull  distribution  (Reference  1), 

Fr(t)  -  PITst]  -  1  — exp[  — (t/d0)"“l 

(83) 

in  which  a0  and  fl0  are,  of  course,  related  to  the  values  of  E[T]  and  <r,r.  Approximate  distribution 
functions  of  T  thus  computed  are  shown  in  Figure  57  as  solid  curves  for  a„  =  0.01,  0.02  and  0.04 
inch,  respectively.  Also  displayed  in  Figure  57  as  circles,  triangles  and  rectangles  are  the 
corresponding  distributions  of  the  test  results  obtained  from  Figure  60.  It  is  observed  from 
Figure  57  that  on  the  basis  of  the  approximation,  Equation  (83),  the  test  results  correlate  very 
well  with  the  present  statistical  fatigue  crack  propagation  model. 


For  additional  comparison,  the  above  procedure  has  been  applied  to  another  set  of  date,  to 
be  referred  to  as  XWPB  where  X  signifies  a  16%  load  transfer  in  the  fasteners  and  WPB  has  the 
same  meaning  as  before.  This  second  data  set  also  has  been  censored  to  include  only  those  with 
a  crack  length  growth  from  0.004  to  0.07  inch,  resulting  in  a  total  of  22  specimens.  The  crack 
propagation  time-histories  of  these  specimens  are  shown  in  Figure  58.  The  linear  regression 
analysis  for  the  estimation  of  model  parameters,  illustrated  in  Figure  59,  resulted  in 

b  -  0.985,  Q0/ab  -  2.4414  x  10  "  \  a^x«>  *  0.12896 
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Figure  57,  Comparison  Between  Weibull-Type  Approximation  for  the  Distributions  of 
Random  Time  to  Reach  a  Given  Crack  Size  and  Actual  Test  Results  for 
WPB  Fastener  Holes 

An  application  of  Equations  (81)  and  (82)  yielded 
M  -  1.0461,  8  -  0.0603 

The  computed  E[T]  ±  <rT  curves  for  a  number  of  A  values  are  shown  in  Figure  60.  When  they 
are  compared  with  the  specimen  mean  and  specimen  standard  deviation,  shown  in  Figure  61,  the 
best  agreement  is  obtained  for  A  «  8000.  Approximate  distributions  of  the  random  times  to 
a0  =*  0.008, 0.026  and  0.07,  based  on  a  Weibull  form,  Equation  (83),  are  shown  in  Figure  62,  along 
with  the  sample  distributions.  Computer  generated  simulations  using  A  -  8000  and  X  =  0.1  are 
shown  in  Figure  63.  Again,  excellent  theoretical  and  experimental  correlations  are  seen  in  this 
Becond  example. 
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Figure  60.  Theoretical  Mean  and  Standard  Deviation  of  Random  Time  to  Reach 

Various  Crack  Sizes  Computed  for  Different  A  Values  for  XWPB  Fastener 
Holes 


Figure  61.  Mean  and  Standard  Deviation  of  Random  Time  to  Reach  Various  Crack 
Sizes  Computed  Directly  from  Some  Actual  Time-Histories  of  XWPB 
Fastener  Holes 
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Figure  62.  Comparison  Between  Weibull-Type  Approximation  for  the  Distribution  of 

Random  Time  to  Reach  a  Given  Crack  Size  and  Actual  Sample  Distribution 
of  XWPB  Fastener  Holes 


Figure  63.  Simulated  Sample  Functions  of  Crack  Propcgation  Time-History  for  XWPB 
Fastener  Holes 
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SECTION  VII 


CONCLUSIONS 


1.  Spectrum  Loading 

A  fracture  mechanics-based  statistical  model  developed  in  Reference  1  for  the  prediction 
of  fatigue  crack  growth  behavior  under  any  single  test  condition  has  been  extended  to  the  case  of 
spectrum  loading.  The  distributions  of  propagation  life  to  reach  any  given  crack  size  and  the 
crack  size  at  any  specific  number  of  load  cycles  were  obtained.  A  specific  block  loading  was 
chosen  and.  applied  to  IN100  laboratory  test  specimens  to  obtain  statistically  meaningful  crack 
propagation  test  data  for  the  verification  of  the  statistical  model,  The  same  block  loading  was 
used  for  the  theoretical  prediction  of  the  crack  propagation  behavior  of  IN100  baaed  on  the 
model.  Comparison  of  the  theoretical  results  with  the  verification  test  data  indicated  a 
reasonable  correlation.  The  model  provided  conservative  predictions  of  both  the  mean  behavior 
and  the  variability  of  the  data.  Thus,  the  proposed  fracture  mechanics-based  statistical  model  is 
quite  practical,  since  it  requires  only  a  small  nonhomogeneous  data  base  for  predicting  crack 
propagation  under  spectrum  loading.  Crack  growth  retardation  or  acceleration  effects  due  to  the 
applied  spectrum  (temperature,  loading)  must  be  addressed  in  further  work  on  this  model. 
However,  the  theoretical  model  is  judged  appropriate  for  practical  applications  due  to  its 
simplicity  and  conservative  predictions. 

2.  Related  Developments 

A  statistical  modol  based  on  the  Paris  crack  growth  rate  function  for  fatigue  crack  growth 
in  IN100  was  investigated  and  applied  to  test  results  at  various  temperatures,  loading 
frequencies,  and  stress  ratios.  Again,  the  distributions  of  life  to  reach  any  specific  crack  size  and 
the  crack  size  at  any  specific  number  of  cycles  were  obtained.  Homogeneous  test  environments 
were  assumed  to  avoid  excessive  extrapolation  into  the  region  where  crack  growth  data  did  not 
exist.  The  life  integration  was  then  based  on  the  best-fitted  crack  growth  rate  parameters.  This 
approach  was  necessary  in  order  to  obtain  the  homogeneous  data  sets,  referred  to  as  the 
extrapolated  teat  results,  for  the  correlation  study. 

The  correlation  between  the  statistical  model  and  the  extrapolated  test  results  was 
reasonable.  Care  should  be  taken  in  the  application  of  this  model  as  the  Paris  crack  growth  rate 
function  is  applicable  only  to  a  certain  region  of  the  stress  intensity  range.  This  model  is 
mathematically  simple  and  practical  for  engineering  applications. 

Closed  form  solutions  have  been  obtained  for  the  statistical  moments  of  the  random  time 
when  a  dominant  crack  reaches  a  given  size  for  a  rather  general  class  of  crack  growth 
mechanisms.  The  key  to  this  success  was  the  approximation  of  the  random  crack  size  by  a 
diffusive  Markov  process.  Theoretically,  the  approximation  introduced  an  error  associated  with 
the  possibility  for  the  crack  propagation  rate  to  be  negative  at  times.  However,  in  terms  of 
statistical  properties,  the  error  was  negligible  if  the  tendency  for  drift  dominates  tho  tendency 
for  diffusion.  The  application  of  the  proposed  theory  to  two  real  examples  seemed  to 
substantiate  our  contention.  The  theoretical  results  have  been  shown  to  correlate  well  with  the 
experimental  results  when  the  parameters  for  the  mathematical  model  were  obtained  from  u 
linear  regression  procedure.  Since  the  calculation  of  a  few  key  parameters  using  such  a 
procedure  doss  not  require  a  very  large  data  base,  the  proposed  theory  is  quite  practical  in  view 
of  the  limited  experimental  results  available  at  the  present  time. 
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36  111000.  29.91  26.02  8.293E-04 

37  112000.  30.15  26.64  2.476E-04 

38  113000.  31.84  27.83  1.689E-03 

39  114000.  38.34  34.10  6.500E-03 

40  115000.  40.30  44.00  1.962E-03 
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